.<0 J- 


(NASA-CK- 1 ?; 7 9 C - V v) 1 - ') Ms M.YUCA 1 K, 0 E 1 1 1. 0 
CF OF: r A 7 1 N s'. 0 fc A A CT c 1 1 3 1 1 C £ r l 

f :<{■: k : :< : :•» o- r f v ? :• si a ; cm.-?ir fin:. g 

tA..iF :G £i. VOl.iJM- I: ^iiALYElH ANC IMILIS 'Ji.CtdS 

lilid 1 Aei-oit ( Ui: LtCsl T< chuclC ‘1 ir £ Nscdi:^ a G J/j** w I Jj4 

v» 

NASA CR-1 67990 
R82-91 5362-40 


Analytical Modeling of Operating 
Characteristics of Premixing-Prevaporizing 
Fuel-Air Mixing Passages 

Vol. I Analysis and Results 


^Olof L. Anderson 
Louis M. Chiappetta 
David E. Edwards 
John B. McVey 



February 1982 


Prepared for 

NASA Lewis Research Center 
Under Contract NAS3-21269 


V 




i* 


- fct- 

. nopon pw, 


4. Tltlt and Subtftla 



Analytical Modeling of Operating Characteristics of 
Premixing-Prevaporizing Fuel-Air Mixing Passages 


7. Author!.) 


0. L. Anderson, L. M. Chiappetta 
D . E . Edwards , J . B . McVey 


I Deport Out 
February 1982 




Vi rWlWrmP^I Uv^WmW ww$mml *w, 

R82-915362-40 


tO. Work Unit No. 


11. Contract or Grant No. 

NAS3-21269 


13. Typo of Raport and ftriod Cov«r*H 
Final Report (Volume I) 

-User's Manual (Volume II ) 

14. Sponsoring Agency Coda 


9. Performing Organization Nama and Addrasi 

United Technologies Research Center 
Silver Lane 

East Hartford, CT 06108 


12. Sponsoring Agency Nama and Address 

NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland, OH 44135 


15. Supplementary Notes 

NASA Project Manager - C. J. Marek 
UTRC Project Manager - 0. L. Anderson 


16. Abstract 

A model for predicting the distribution of liquid fuel droplets and fuel vapor in pre- 
mixing-prevaporizing fuel-air mixing passages of the direct injection type is reported 
herein. This model consists of three computer programs; a calculation of the two- 
dimensional or axisymmetric air flow field neglecting the effects of fuel a calcula- 
tion of the three-dimensional fuel droplet trajectories and evaporation rates in a 
known, moving air flow; a calculation of fuel vapor diffusing into a moving three- 
dimensional air flow with source terms dependent on the droplet evaporation rates. The 
fuel droplets are treated as individual particle classes each satisfying Newton's law, 
a heat transfer, and a mass transfer equation. This fuel droplet model treats multi- 
component fules and incorporates the physics required for the treatment of elastic 
droplet collisions, droplet shattering, droplet coalescence and droplet wall interac- 
tions. The vapor diffusion calculation treats three dimensional, gas-phase, turbulent 
diffusion processes. The analysis includes a model for the autoignition of the fuel- 
air mixture based upon the rate of formation of an important intermediate chemical 
species during the pre-ignition period. This species is produced both within the 
vicinity of the fuel droplets and throughout the diffusing fuel vapor-air mixture. 

The model, as represented by these computer codes, is applied to two premixing fuel- 
air mixing passage designs and the results are discussed. An application of the 
autoignition model is also presented. 


17 Key Word! (Suggested by Author(s)) 
Premixing-Prevaporiz ing Passages 
Autoignition Model 
Two-phase Flow 

Two- and Three-Dimensional Flow Models 
Computer Program 


19 Security Qe»if (of tM report) 20 Security CNoif (of tW» 

Unclassified Unclassified 


18. Distribution Statement 


21. No. of 



* For sale by the National Technical Information Service, Springfield. Virginia 22151 
















R82-915362-40 


Analytical Modeling of Operating Characteristics of 
Premixing-Prevaporizing Fuel-Air Mixing Passages 


Volume I - Analysis and Results 


TABLE OF CONTENTS 


1.0 SUMMARY. 


2.0 INTRODUCTION 


3.0 ANALYSIS OF AXISYMMETRIC AIRFLOW 


3.1 General Approach. . 

3.2 Coordinate System . 

3.3 Equations of Motion 

3.4 Turbulence Model. . 

3.5 Numerical Methods . 

3.6 References 

3.7 List of Symbols . . 


4.0 ANALYSIS OF FUEL DROPLET BEHAVIOR . 


4.1 General Approach 

4.2 Momentum Equations 

4.3 Heat and Mass Transfer Equations 

4.4 Evaluation of Thermodynamic and Transport Properties 

4.5 Evaluation of Distillate Fuel Effects 

4.6 High Pressure Effects 

4.7 Droplet Class and Distribution Functions 

4.8 Droplet Collision Model 

4.9 Droplet Shattering Model 

4.10 Droplet-Solid Boundary Interaction Model 

4.11 Numerical Methods 

4.12 References 

4.13 List of Symbols 

4.14 Tables 

4.15 Figures 


1 


i 


R82-915362-40 


TABLE OF CONTENTS (Cont'd) 


5.0 ANALYSIS OF DIFFUSION OF FUEL VAPOR. 


General Approach 

Fuel Vapor Diffusion Equation 
Numerical Methods ...... 

References 

List of Symbols 


6.0 COMPUTER MODEL CALIBRATION 


Calibration of Fuel Droplet Model . . 

Model Sensitivity Study 

Calibration of Vapor Diffusion Model. 

References 

Tables 

Figures 


7.0 PARAMETRIC STUDIES 


Swirl Tube Premixing Passage. . 
Series Staged Premixing Passage 

Reference 

List of Symbols 

Figures 


8.0 ANALYSIS OF AUTOIGNITION 


8.1 General Approach 

8.2 Basic Equations 

8.3 Autoignition Models 

8.4 Numerical Methods 

8.5 Results and Discussion 

8.6 References 

8.7 List of Symbols 

9.0 CONCLUDING REMARKS 
APPENDIX A - LITERATURE SURVEY 


11 


R82-915362-40 


% 


V 


Analytical Modeling of Operating Characteristics of 
Premixing-Prevaporizing Fuel-Air Mixing Passages 

Volume I - Analysis and Results 


1.0 SUMMARY 


A model for predicting the distribution of liquid fuel droplets and fuel 
vapor in premixing-prevaporizing fuel-air mixing passages of the direct injection 
type is reported herein. This model consists of three computer programs: a 

calculation of the two-dimensional or axisymmetric air flow field neglecting the 
effects of fuel; a calculation of the three-dimensional fuel droplet trajectories 
and evaporation rates in a known, moving air flow; and a calculation of fuel 
vapor diffusing into a moving three-dimensional air flow with source terms depen- 
dent on the droplet evaporation rates. The air flow calculation can treat com- 
pressible swirling flows in arbitrary ducts with arbitrary pressures, temperatures 
and velocities as initial conditions. The fuel droplets are treated as individual 
particle classes each satisfying Newton’s law, a heat transfer, and a mass trans- 
fer equation. Each particle class has a number density such that summation over 
all particle classes yields the fuel flow rate. This fuel droplet model treats 
multicomponent fuels and incorporates the physics required for the treatment of 
elastic droplet collisions, droplet shattering, droplet coalescence and droplet 
wall interactions. The vapor diffusion calculation treats three-dimensional, gas- 
phase, turbulent diffusion processes with the turbulence level determined by the air 
flow calculation and the source terms determined by the droplet evaporation rates. 

The analysis includes a model for the autoignition of the fuel-air mixture 
based upon the rate of formation of an important intermediate chemical species 
during the pre-ignition period. This species is produced both within the vicinity 
of the fuel droplets and throughout the diffusing fuel vapor-air mixture. Since 
chemical reaction rates may depend upon the local mixture temperature, the local 
mixture temperature is adjusted for the effect of fuel evaporation. 

A preliminary calibration of the computer codes and a parametric study 
showing the effects of the air flow conditions and initial droplet conditions on 
the evaporation rate are included. The model, as represented by these computer 
codes, is applied to two premixing fuel-air mixing passage designs and the re- 
sults are discussed. An application of the autoignition model is also presented. 
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2.0 INTRODUCTION 


One combustion control strategy for meeting Governmental regulations on 
emissions of pollutants from internal combustion engines which has recently 
received considerable attention is the use of the premixing, prevaporizing combus- 
tion concept whereby uniform, homogeneous fuel-air mixtures are delivered to the 
combustion chamber in such proportions that the gas temperature-time history per- 
mits complete oxidation of the hydrocarbon fuel but does not permit significant 
production of the oxides of nitrogen. Methods of achieving premixed, prevaporized 
fuel-air mixtures include external vaporization schemes whereby the fuel is vaporized 
before being mixed with air and direct injection of a finely atomized spray into 
the airstream. This effort is concerned with the prediction of the distribution 
of liquid and vapor fuel produced by the direct injection method. 

Distribution of the liquid fuel throughout an airstream is primarily effected 
through the use of multi-point fuel injection, by imparting velocity components to 
the fuel droplets which are normal to the airstream velocity components, and by 
turbulent diffusion of the fuel throughout the airstream. Vaporization of the 
spray is achieved by heat transfer from the airstream to the droplets and mass 
transfer of fuel vapor from the droplet surface to the droplet environment. Knowl- 
edge of the magnitude of the vaporization time is of particular concern in the case 
of gas turbine engines wherein fuel is injected into high temperature compressor 
discharge air due to the requirements that the vaporization time be significantly 
less than the autoignition time of the fuel-air mixture and that the required passage 
length be consistent with engine dimension limits. From a design point of view, 
what is required is an analytical procedure whereby the atomization and initial 
distribution characteristics which a given fuel injection system will deliver can 
be evaluated in terms of the fuel distribution and degree of vaporization which can 
be achieved within a given mixer/vaporizer section length. 


For prevaporizing premixing fuel-air passages, autoignition must be prevented 
because the uncooled engine hardware in this duct can undergo catastrophic fail- 
ure as a result of the sudden large heat release following autoignition. Clearly, 
exploitation of the concept requires that the residence time within the fuel pre- 
paration passage be long enough to achieve essentially complete vaporization and 
yet short enough to preclude the occurrence of autoignition. An analysis is re- 
quired, therefore, to predict if autoignition occurs in premixing passages at 
simulated gas turbine engine operating conditions. 



This report describes the analytical models which were developed to predict the 
operating characteristics of premixing-prevaporizing fuel-air mixing passages of 
the direct injection type. The technical approach adopted for developing these 
models is to separate the problem into three parts each with its own computer code. 
These three parts are: calculation of the two-dimensional or axisymmetric gas flow 

field (ADD code), calculation of the three dimensional nonequilibrium heating and 
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vaporization of the fuel droplets (PTRAK code), and calculation of the three-dimensional 
turbulent diffusion of fuel vapor in air (VAPDIF code). In applying this approach, 
it is implicitly assumed that the air flow behavior can affect the fuel droplet 
behavior, but the fuel droplet behavior does not affect the air flow behavior. This 
is justified because the mass fraction of fuel droplets and fuel vapor is small. 
Similarly, it is assumed that the fuel droplet behavior affects the fuel vapor 
behavior but the fuel vapor behavior does not affect the fuel droplet behavior. This 
decoupling (or weak interaction) assumption of the problem allows a much simpler 
solution wherein the air flow behavior, fuel droplet behavior, and vapor diffusion 
can be analyzed in succession. 

This report also describes a model for predicting autoignition in two-phase 
turbulent flows. It is assumed in the model that autoignition is determined by 
the rate of formation of an intermediate chemical species, it is also assumed 
that production of this species has a negligible effect on fuel vapor concentra- 
tion. Sources of this species are calculated by the PTRAK code within the vicinity 
of the fuel droplets; additional sources are calculated by the VAPDIF computer 
program as the fuel vapor diffuses throughout the premixing passage. These sources 
are then used in the VAPDIF code to determine the local concentration of the rate- 
controlling species. A criterion is applied to determine if autoignition of the 
fuel-air mixture has occurred. 


Some chemical reaction rates are functions of the local gas temperature. A 
first-order correction is made to the temperature distribution due to evaporation 
of the fuel prior to determination of the local rate of production of the inter 
mediate species. Consistent with the weak interaction assumption described 
earlier, the effect of this temperature correction on the flow field or vaporiza- 
tion rate is assumed to be negligible. 

The model, composed of the three computer programs, 1) treats the three dimensional 
behavior of the fuel droplets and fuel vapor in a moving gas stream, 2) represents a 
reasonable compromise between rigor and empiricism, and 3) yields practical results 
with a reasonable expenditure of computer time. This report describes the analytical 
models incorporated into the computer codes and presents the results of a preliminary 
calibration of the models, a sensitivity study, and an analysis of two premixing 
passage designs that illustrate the features of the analytical models incorporated 
into the computer programs. An application of the autoignition model is also pre- 

rented . 
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3.0 ANALYSIS OF AXISYIMRkIC AIRFLOW 


3.1 General Approach 

The ADD code was developed to solve the internal flow weak interaction problem 
using a forward inarching numerical procedure that does not require an interaction 
between the inviscid core flow described by the elliptic Euler equations and the wall 
boundary layers described by the parabolic boundary layer equations. The basic 
mathematical description and justification of the technical approach is given by 
Anderson (Ref. 3.1). A more detailed derivation and description of the various 
features of the computer code is given in Ref. 3.2. Applications of the computer 
code are given in Refs. 3.3 and 3.4. 

This method can be described in the following manner. First, an orthogonal 
coordinate system is constructed for the duct from the potential ^owsolutionsuch 
that the stream function forms the coordinate normal to the wall and the velocity 
potential forms the coordinate tangent to the wall. Since the potential ow 
streamlines approximate the real streamlines, the equations of motion may be 
greatly simplified by assuming that the velocity normal to the potential flow 
streamlines is small compared to the streamwise velocity. This P*® 6 **"® *****?* 
the governing viscous flow equations to a parabolic system of partial differentia 
equations which can be solved by a forward marching numerical integration procedure. 
Furthermore, it can be shown (Ref. 3.1) that the resulting solution has the same 
order of accuracy for viscous flows as the streamline curvature method has for 
inviscid flows. Thus, the inviscid-flow weak interaction with the waU boundary 
layer is solved without the need for iterations between different flow fields, 
i.e., the inviscid core flow and the boundary layer flow. 

Symbols used in this section are defined on pages 17-18. 
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3.2 Coordinate System 

The equations of motion are solved in a streamline orthogonal coordinate system 
(n, s, <t>). This coordinate system is generated by first solving for the plane 
potential flow through a duct with the same cross-section that the annular duct 
makes with the meridional plane. The normal coordinate n is the stream function and 
the streanwise coordinate s is the velocity potential. Rotation about the axis of 
symmetry produces an axisymmetric orthogonal coordinate system uniquely suited to 
solve the problem. 

The (n, s) coordinates are related to the physical coordinates (r, z) through 
LaPlace' s equation. 


dz 2 


(3.2.1) 


(3.2.2) 


The metric scale coefficients are the same in both directions and are equal to the 
inverse of the magnitude of the potential flow velocity. 




(3.2.3) 


Lengths along the streamlines and potential lines are given by 


dx = ds/v 


(3.2.4) 


dy * dn/v 


(3.2.5) 


and curvatures of streamlines and potential lines are given by 
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(3.2.6) 

(3.2.7) 

If the duct wall contours r w (z) are specified, Eqs. (3.2.1) and (3.2.2) are solved 
using the complex potential and the Schwartz-Christoffel transformation (Ref. 3.5) 
using the method described in Ref. 3.1). 


r t dn 

JL.» 

r n d* 
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3.3 Equations of Motion 


The ADD code equations* given below* are written in an orthogonal streamline 
coordinate system where n is the normal coordinate (potential flow stream function) 
and s is the streamwise coordinate (potential flow velocity potential) . The metric 
scale coefficient is the same in both the n and s directions and is equal to (1/V) 
where V is the magnitude of the potential flow velocity. 


d* r*U # 


an 


(3.3.1) 
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(3.3.7) 
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(3.3.8) 


P * pRT 


(3.3.9) 


1-4 • Cp ln(T/T 0 ) - Rin (P/P.) 


(3.3.10) 
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Equations (3.3.1) through (3.3.10) form a sat of sight first order partial 
differential equations and two algabrale equations which may be used to solve for 
ten unknowns. The boundary conditions for this problem are given by 


u t (0,s)=0 


for the ID wall and 


u+(o,s)=o 


QnK), s) = 0 


< ff (0, s) s 0 


u s O,s) = o 


U^(l,s) = 0 


Qn(l,s)-0 




for the OD wall* 
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3 .4 Turbulence Model 

Three turbulence models have been programmed into the ADD code. These 
turbulence models are: (1) eddy viscosity model (Ref. 3.1), (2) modified eddy 

viscosity model (Ref. 3.6) which uses the streamline curvature corrections of Bide 
(Ref. 3.7) and Bradshaw (Ref. 3.8), and (3) two equation turbulence model (Ref. 3.6) 
which is based on the work of Chen (Ref. 3.9) using the curvature corrections of 
Launder, et al- (Ref. 3.10). 


Eddy Viscosity Model 

The eddy viscosity model is essentially an equilibrium model based on Prandtl’s 
mixing length theory (Ref. 3.11). The outer layer and freestream flow have an eddy 
viscosity which varies with s only. This freestream eddy viscosity is given by. 




ft. U m 
2 


h 


•/('-£> 


(3.4.1) 


where x is an empirical constant assumed to be 0.016 (Ref. 3.12) and P^U® is taken 
to be the maximum in the duct. Hence, for thin boundary layers, the integral in 
Eq. (3.4.1) reduces to 2 6* where 5* is the displacement thickness. Alternately 
with large freestream distortion, Eq. (3.4.1) can be thought of as: 


ft T *X p m h (u.-u) 


(3.4.2) 


where U is the average velocity and h is the duct length. Equation (3.4.2) is a 
typical wake mixing length model. 

For the "wall layer", the turbulence model derived by VanDriest (Ref. 3.13) is 
used . 


/) 2 [i-exp(-yVA*)] 


where U + , Y + are universal boundary layer coordinates given by. 




i y / r m f P //i-w 

where the subscript w refers to values at the wall. 
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The genc'.lly accepted values of the empirical constants x» 1 are taken to be 
0.016, 0 . 40 . 26.0, respectively. For swirling flows U and U+ are determined from 
the resultant velocity. 


Modified Eddy Viscosity Model 


The modified eddy viscosity model Is a two layer model where the Prandtl mixing 
length l is defined by 


/V' 1 ’ I* 


(3.4.6) 


In the wall region, the mixing length is given by 

p U*Jt 

“7^- : «Y + [l-6Xp(-Y + /A + )] 


(3.4.7) 


In the outer region or wake region Including the freestream, the mixing length is 
Eq. (3.4.6) 


H- t s P A{[J 00 - u) 


(3.4.8) 


where 


l = Xh 


(3.4.9) 


Eide and Johnston (Ref. 3.7), modified the mixing length in the wall layer by a 
factor F given by 


F = «/V + 0, R i. + V* 


(3.4.10) 


where d s and 6^ are empirical constants and R 1<( and Rj^ are the Richardson numbers 
for streamline curvature and swirl defined by 


; H ( ,= 2S,(I+S,I 

(3.4.11) 

; % !2 V + V 

(3.4.12) 


Elde estimated that 0 8 - 0^ • 6.0. 
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The outer region of the flow field .8 governed by large scale mixing. Using 
the suggestion of Bradshaw (Ref. 3.8) for a bulk Richardson number, the outer layer 
mixing length is modified by a factor F 


7=1 + + a^)(2Y/h-.) 


(3.4.13) 


where 


Sis 1 - h/r s 


(3. A .14) 




h u<£ 
s 


(3.4.15) 


Bradshaw estimates values of a « a. * 6.5. 

s (J> 

Two Equation Turbulence Model 

The two equation turbulence model programmed into the ADD code is based on the 
work of Chen (Ref. 3.9) with the curvature corrections of Launder (Ref. 3.10). In 
the ADD code coordinates, these equations are given by 
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The empirical constants are given by: 


C, = 1.35 
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<r k * 1.0 


<r« =1.3 


Finally, the boundary conditions are given by 


k(0) = «(o) = 0 


kli) * «(i) * 0 
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3.5 Numerical Methods 

With the relationship between |i£ and the naan flow specified, Eq. (3.3.1) 
through (3.3.10) can be solved by a forward marching numerical integration schema. 
Equations (3.3.1 ) through (3.3.10) are first linearised by expanding all dependent 
variables in a Taylor series expansion in the marching direction (s) , and terns of 
0 (As 2 ) are dropped. Finite difference equations are then obtained using the two 
point centered difference scheme of Keller (Refs. 3.1A, 3.15). The resulting 
matrix equations are (10 x 10) block tridiagonal and are solved by block factorisa- 
tion using the method of Varah (Ref. 3.16). 

The numerical solution is second order accurate in the n direction, first 
order in the s direction, and linearly stable. The As step sise is limited not by 
linear stability conditions but by the required accuracy in the Taylor series 
expansion in s. 
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3.7 List of Symbols 

Van Driest constant (26.0) 

Specific heat of air 

Constant in Two Equation Turbulence Model 
Rate of Strain 

Ratio mixing length to flat plate mixing length 
Ratio free stream mixing length to wake mixing length 
Duct height 
Entropy 

Turbulence kinetic energy 



* j 

Mixing length 

Mixing length flat plate 

] 

Normal coordinate (stream function) j 

Static pressure 

Turbulence production I 

Rad ius j 

Radius of curvature 
Gas constant 

Richardson numbers for streamline curvature 

Turbulent Richardson number \ 

Richardson number for swirl j 

Turbulent Reynolds number 
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List of Symbols (Cont'd) 

Streamwise coordinate (velocity potential) 

Temperature 
Velocity components 
Maximum velocity 
Average velocity 
Universal velocity 
Friction velocity 

Reciprocal of metric coefficient (potential flow velocity) 

Distance along s and n coordinates 

Universal distance from wall 

Axial distance 

Empirical constants 

Empirical constants 

Displacement thickness 

Turbulence dissipation 

Prandtl constant (.41) 

Molecular viscosity 
Effective viscosity 
Turbulent viscosity 
Density 
Stress 

Clauser constant (.016) 

Stream function 
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4.0 ANALYSIS OF FUEL DROPLET BEHAVIOR 


4.1 General Approach 

The behavior of the fuel spray Is calculated using the Particle Tracking (PTRAK) 
computer program which is described in this section. 

The analysis assumes lean equivalence ratios, i.e., that the mass fraction of 
fuel is small. Under these conditions, the air flow behavior can be solved first 
using the analysis described in Section 3.0 and incorporated into the ADD code. 

The pressure, temperature and velocity field are then stored on a data file and 
used in the fuel droplet analysis. Thus the air flow behavior can effect the fuel 
droplet behavior, but the fuel droplet behavior does not effect the air flow 
behavior. 

The spray is divided into a number of classes (up to 1250). In each class, 
droplet diameter, velocity components, and position in the flowfield are specified; 
these parameters may differ among the various classes. The PTRAK model calculates 
the trajectory, vaporization rate, and temperature variation of each class of drop- 
lets. Large droplets may shatter into smaller droplets. Droplets from two classes 
may coalesce into larger entities. Droplets may collide with the duct walls and 
either rebound or undergo additional vaporization due to heat transfer to the drop 
let from the wall. The fuel may be either a pure compound or a distillate liquid. 

As a result of these processes, a distribution of fuel vapor sources is calculated 
for use in the vapor diffusion computer program (see Section 5.0). 

Symbols used in this section are defined on pages 57-60. 
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4.2 Momentum Equations 

The equations of motion for the fuel droplet trajectories presently coded into 
the PTRAK computer program are given below. It should be noted that these equations 
of motion are written in the orthogonal streamline coordinate system used by the 
ADD code. The additional terms on the right hand side are acceleration terms 
generated by the curvilinear coordinate system. V g , V n , and are the fuel drop- 
let velocities and s, n , <p «e the coordinates of the fuel droplet in the ADD code 
coordinate system. The drag of the fuel droplet is written in terms of a drag 
coefficient, Cp. 



PCq A 

<Jt * 


jL , v„v t 


dv 

ds 


(4.2.1) 


(4.2.2) 


l.ii 


4i.*vv t 


dt 

dn . 
dt 


W r 


*L t vt 

dt t 

where the cross-sectional area of the fuel droplet is given by 


(4.2.4) 

(4.2.5) 

(4.2.6) 


A = irD /4 


(4.2.7) 


The drag coefficient is determined from correlations reported by Dickerson and 
Schuman (Ref. 4.1): 


C D = 27 Re 


* .84 


Co s 0.271 Reo 
C 0 - 2 0 


0217 


0 < Reo <80. 
80 < Reo < I0 4 


(4.2.8) 

where the droplet Reynolds number is based on film properties and the relative velocity. 


Re 0 > 10' 


"«0 ~ Pm DAU/Mm 


(4.2.9) 
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4.3 Heat anc 1 Mass Transfer Equations 

In the case of a (single component) droplet vaporizing in a gas stream of 
uniform temperature, the droplet temperature history can be described in terms of 
a heat-up period, during which time the droplet temperature changes, and an 
"equilibrium vaporization" period, during which time the droplet temperature 
remains constant. The physical processes of major significance during these 
periods are the transport of heat to the droplet and the transport of vapor away 
from the droplet surface. The net difference between the energy flux to and 
away from the droplet accounts for changes in the droplet temperature. 

The rigorous mathematical treatment of the droplet vaporization and heating 
problem requires the solution of time-dependent partial differential equations for 
the diffusion of mass and energy for conditions both within the droplet and in the 
surrounding gas-phase. Understandably, a great amount of insight into the behavior 
of a droplet can be obtained from such an approach. However, there are several 
reasons why such a procedure is not incorporated into the PTRAK code. First, 
solutions to the set of partial differential equations would have to be obtained 
for hundreds of droplet classes throughout the premixing passage; computer run 
times would become unacceptably large. Second, it is not clear whether sufficient 
information is available for calculating heat and mass transfer coefficients for 
droplets in forced convection (due to the relative velocity between the droplet 
and the air). Third, there is a paucity of data useful for verifying a model of 
this complexity especially if the model includes the effects of distillate fuels on 
heat and mass transfer rates. 

Instead of using a rigorous mathematical treatment of the droplet heat and 
mass transfer processes in the development of the PTRAK program, an existing semi- 
empirical model was extended for use with distillate fuels. The model described 
herein is based largely on the work of El Wakil, et.al (Ref. 4.2) and Priem and 
Heidmann (Ref. 4.3). 

There are three assumptions that make the mathematical treatment of droplet 
heat and mass transfer quite tractable. The first is the quasi-steady assumption 
by which the transport processes are assumed to be steady at every instant; that 
is, at each time t during droplet vaporization and heating, the fluid and thermo- 
dynamic conditions of the droplet and gas system adjust instantaneously to steady- 
state conditions. The quasi-steady assumption is not an equilibrium assumption. 
With a quasi-steady assumption, the droplet temperature and vaporization rate 
may change with time. With an equilibrium assumption, the droplet temperature is 
constant (and equal to its wet bulb value); an equilibrium (constant) vaporization 
rate does not exist but it can be shown that the ratio of vaporization rate to 
droplet diameter is nearly constant at equilibrium. 

The second assumption is that conditions within the droplet are uniform. This 
permits the governing partial differential equations for determining droplet 
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temperature and mass (or diameter) to be reduced to extremely simple forms for 
quick solution. 

The third assumption is that both the droplet and its surroundings are spheri- 
cally symmetric. Obviously, the relative velocity between the droplet and the gas 
causes the droplet to change shape and alter the characteristics of the velocity, 
thermal, and mass diffusion boundary layers. Correlations developed for vaporizing 
droplets in a forced convection environment provide the appropriate transport 
coefficients. 

Using these assumptions, the governing partial differential equations for the 
gas phase are reduced to a set of steady-state, ordinary differential equations with 
known solutions. These solutions, together with the simple equations for the liquid 
phase, are combined in a forward-marching time integration procedure (the time- 
equivalent of the spatial integration of the droplet equations of motion. Sec. 4.2). 

The geometry of the system under consideration is shown in Fig. 4-1. The 
droplet is surrounded by a boundary layer of uniform thickness and properties. At 
the outer edge of this boundary layer, gas conditions correspond to those of the 
flow field calculated by the ADD code and will be referred to as the "air" conditions. 
The convention adopted in Ref. 4.2 and retained here is that the thickness of the 
boundary layer is equal to the instantaneous droplet radius. For a complete develop- 
ment of the equations for vaporization and heating rate, the reader is referred to 
Ref. 4.2 and 4.3. The vaporization rate expression has been modified (Ref. 4.4) 
to account for a finite amount of fuel vapor in the gas outside of the boundary 
layer. The effect of this "back diffusion" tern is negligible for lean systems; 
it has been retained for use in the PTRAK code but is presently rendered inoperable. 

The equation for the vaporization rate is: 


*1 = A S K Pf,s a 


(4.3.1) 


where 


0 ’ , £«_ |n 

P '.» \ P o D ‘,s / 


(4.3.2) 


The term p f , w is the partial pressure of fuel vapor outside of the boundary layer. 
(„ = o in* the present version of PTRAK.) The determination of the mass transfer 

coefficient is presented in the next section. The droplet heating rate equation is: 


dT L _ 

dt 


ds - W L X 

■f ° 3 Pl C PL 


(4.3.3) 


23 


R82-915362-40 


ORIGINS pA ®f 
OF POOR QUALITY 



where the total amount of heat arriving at the droplet surface la: 


q s * hA s ( T 0 -Tt) e * 2 _, 


and the parameter z is: 


(4.3.4) 


*1 

hAs 


(4.3.5) 


The droplet surface area, A 8 , is calculated from: 


A« = fD 


(4.3.6) 


,l,o de terminal ton of tho heat transfer coefficient. h, la presented la the next 
section. Equation (4.J.1) states that some of the heat arriving at the dr p 
surface Is used to vapor, to an amount of fuel, i,, while the remainder Is used 
TC at the droplet. At cgulllhrlum. g a la e,ua, to w » and thedropet 
(the wet-bulb temperature) is constant (for a pure fluid). Equation (4.3.4) is 
product of the forced-convection heat transfer rate for a nonvaporizing system, 
hAs (To - T L ), and a correction term for the energy carried away from dm 

to mass transfer (t.e.. "bloving" of the boundary layer). t/(e - D • ' Thin term 
approaches unity for a tero mas. transfer rate and tern for an Infinite mas. transfer 
rate (no energy reaching the droplet surface) 


At any Instant, the droplet mass is given by 


m L = 


T °V 


(4.3.7) 


Then, since 


dm L 

~3T~ 


- w, 


(4.3.8) 


and assuming that (T,.3 only. Equations (4.3.3). (4.3.7) and (4.3.8) and 

bv l ombliu'ii to y WIs 

2 |. .0* dg L dft| 

(4.3.9) 


<*2. = 

dt *0 Vl 


"»■ * 6 dT. dt 


where d»<j /dTj Is evaluated as described In the next section. 
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4.4 Evaluation of Thermodynamic and Transport Properties 

In this section, the procedures used to calculate the heat and mass transfer 
coefficients are presented together with the techniques used to calculate the 
required thermal properties. The discussion contained herein is limited to flow 
situations in which the local static pressure is well below the critical pressure 
of either the fuel or the air. Most of these procedures are also used in the; 
high pressure model for droplet heating and vaporization. The discussion of the 
effect of high pressure is deferred until Section 4.6. 

(1) Mass Transfer Coefficient: 

The mass transfer coefficient, K, is obtained from a correlation of the form 
developed by Eanz and Marshall (Ref. 4.5) for droplet vaporization: 


= 


KDTmRo 
•^m 771m 


1/2 1/3 

a + b Re D Sc 


(4.4.1) 


The first term represents the mass transfer due to free convection while the second 
term represents the mass transfer due to forced convection. Experimental data are 
required to calibrate suitably the coefficients (a,b). No useful data on the 
vaporization rates of fuel sprays are currently available; therefore the only recourse 
at present is to use the Ranz and Marshall results which were obtained for the 
vaporization of water droplets: 

a = 2.0 
b =■ 0.6 

These coefficients have been incorporated into the PTRAK computer program. 

(2) Heat Transfer Coefficient: 

A companion expression was developed by Ranz and Marshall for the heat transfer 
coe f f ic ient , h : 


m - nD ' „ I/* 1/3 

Nu h " T — = 0 ♦ b R« 0 P r 
* m 


(4.4.2) 


For reasons similar to those presented during the discussion of K, the coefficients 
have been assigned the values; 


a' 2.0 
b' * 0.6 
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in the PTRAK code. 

(3) Dimensionless Transport Numbers; 

The Reynolds, Schmidt, and Prandtl numbers are defined as: 

R «0 = /° m DAU ^m 
Sc = 



Pr - f» m Cpm/ km 


(4.4.5) 


The thermal coefficients used in these parameters are evaluated at a mean boundary 
layer film condition. 

(4) Mean Film Conditions 


At any instant during the evaluation of the vaporization and heating rates, it 
is assumed that the appropriate properties are constant throughout the film 
surrounding the droplet (Fig. 4.1). In this manner, the differential equations 
for heat and mass transfer can be integrated to yield Eq. (4.3.1) and (4.3.4). 

The mean temperature, Tm,, used in the PTRAK code is the log mean value (as 
suggested in Ref. 4.2). 


T T 0- Tl - 

Tm " ln(To/T L ) 


(4.4.6) 


In the development of the droplet vaporization and heating models used in the PTRAK 
program, it was assumed that the air temperature varied from 400K to 900K. If it 
is also assumed that the fuel temperature is typically 300K or greater, then T a / 

T L <_ 3. As a result, the log mean temperature will differ by no more than ten percent 
from the arithmetic mean temperature, 1/2 (T a + T L ) . (If the model were used with 
much higher gas temperatures, such as a fuel droplet might encounter in a combustion 
chamber, then the difference between alternative definitions of mean temperature 
can be quite substantial.) 

It is also assumed that the fuel vapor and the air thermal coefficients are 
obtained at a mean fuel vapor concentration (mole fraction): 


Ym 2 *5“ (Yf f08 4- Yf f ») 


(4.4.7) 
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The mole fraction is also equal to the ratio of partial pressure to local static 
pressure so that: 


Y f,» = Pf,oo /Po 
Y f,s = Pf,s 7 Po 


(4.4.8) 

(4.4.9) 


As noted earlier, it is assumed in the present version of the PTRAK code that Pf,^ 
is zero. The partial pressure of fuel vapor at the droplet surface, Pf» s » is assumed 
to be equal to the vapor pressure of the fuel at the droplet temperature, T . 


The mean molecular weight of the mixture at the mean condition is: 

= Y m ^ f Ml-Y m )7n Q 


(4.4.10) 


The mean mass fraction of the fuel vapor is 


c m = Y m^f /Vf t 


(4.4.11) 


The five mean thermal properties that must be evaluated are density 
(p m ), molecular viscosity (u m ), thermal conductivity (k^,) , heat capacity (Cp m ) , 
and diffusion coefficient G£ m ) . It is assumed that no pressure gradient exists 
across the boundary layer. Then from the ideal gas equation of state: 


Pm ~ Po w lm^ tR o T m* 


(4.4.12) 


(see also the discussion in Section 4.6). The mean molecular viscosity is evaluated 
as a mole weighted average: 


Mm s Y m Mf + C-Y m )^ < 


(4.4.13) 


where the viscosities, v a and y f are evaluated at T m . (A mole weighted average is 
used because viscosity is a manifestation of momentum transfer due to molecular 
collisions and the frequency of collisions is proportional to molecular concen- 
tration). Similarly, since heat transfer by conduction is a function of the 
frequency of molecular collisions, the mean thermal conductivity is: 


27 


R82-91 5362-40 




omsmu- s 

OF POOR QUALITY 

* T m k, Ml-Y m )k 0 


<4.4.14) 


where the thermal conductivities are evaluated at T . 

n 

The heat capacity is a measure of the ability of a system to store energy so 
that the mean heat capacity (following Ref. 4.2) Is: 

Cpm = CmCpf + Cpa (4.4.15) 

where the heat capacities (defined in terms of energy /mass /degree) are evaluated 
at T m . 

The major components of the gas in the film surrounding the droplet are those 
of air under flow conditions likely to be found in premixing passages. The use 
of Eqs. (4.4.13), (4.4.14) and (4.4.15) in lieu of more exact expressions for 
mixtures of fuel vapor and air would have only a small effect on the calculated 
heat and mass transfer rates. The diffusion coefficient, cannot be estimated 

as simply, however, with the same level of precision. As suggested in Ref. 4.2 the 
diffusion coefficient is obtained using the kinetic theory result (Ref. 4.6): 


<0m = 


0.002628 


3 'll If 


~ 2fr f 

Po v m 


(4.4.16) 


where 


<r m = 1/2 ( <t 0 f <r f ) 


ft = 11 ( T ) 


T*= T m /(«/k). 


(4.4.17) 

(4.4.18) 

(4.4.19) 


/JL\ = (—) 

(-L) 

(4.4.20) 


U; m vu/ f 

\ k t Q 


« 


The collision Integral function, 0 , is tabulated in Table 4.1. The force constants 
for air are: 

- 3.689 X 
( k )o * ® 4K 
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Fore* constants are not available for typical jet fuels. It Is reconeadad that 
the force constants for the nolecule ( (C^I^q) (the heaviest paraffin reported in 
Ref. 4.6) be used: o- 


(*), 

(5) Fuel Properties Required in the Model 


8.448 
240 K 


Fourteen properties for each fuel must be determined for use in the PTRAK 
model in the most general problem. The equations described in the preceeding dis- 
cussion account for most of the required properties. The remaining properties are 
needed for special features of the model that will be described below. 

Five of the properties are constants. These are the fuel molecular weight (ty f ) , 
critical pressure (P C f) » critical temperature (T c -) and the force constants (o f ) 
(</k f ). The force constants are generally not available for jet fuels; at present, 
the constants presented above are recommended. 

There are four vapor phase properties. Three of these are input to the 
model in the form of polyn' ial expressions: 


4 N *1 

Property = £ A N T m 
NX 


(4.4.21) 


where "property" is the fuel vapor molecular viscosity (pf ) , the thermal conductivity 
(k f ) or the heat capacity (Cp f ) . The fourth parameter is the vapor pressure (p f>8 ) 
which must be evaluated using a procedure described in Sec. 4.5; it will be shown 
therein that the vapor pressure will be calculated using an expression of the form: 


ln x,»* ♦•ff 


(4.4.22) 


where the "constants" (a n , 6 n ) are determined in accordance with Sec. 4.5. 

Five properties of the liquid phase must also be calculated. The liquid density 

( p .) and heat capacity (c p ) are calculated from: 

L r L 

4 

Pi <<* Cp L > * I A N T"" (4.4.23) 

N*l 


•:; A 


-V:- ■/* 


The liquid molecular viscosity (p^) d«i erm l ne< l from: 


In Pi * I A w T 

N>l 


N-l 


(4.4.24) 


29 


R82-915362-40 


ORIGINAL PAGE IS 
OF POOR QUALITY 


The semi-logarithmic form ie used beceuae the liquid vlacoalty la a strong 
function of temperature. Often, an alternative expression of the fora 


ln H£r* c )] 


a ln*r L ♦ b 


(4.4.25) 


is used for calculating viscosity (Ref. 4.7). The form of Eq. (4.4.24) was selected 
for use in the PTRAK code since it represents the data adeauately and its coefficients 
are determined more conveniently than those used in Eq. (4.4.25). The correlation 
for surface tension (S^) reported in Ref. 4.7 is used: 




= 2.12 ( Tc f - 6 - T L ) 


(4.4.26) 


In this expression, the liquid density is evaluted at T^. The heat of vaporization 
(A) is calculated using the Clausius-Clapeyron equation (e.g.. Ref. 4.8, p. 104) 
which relates A to the slope of the vapor pressure curve: 


dl "Pf,s 

dT 


-xtyf 
R o T l 2 


(4.4.27) 


Then, using Eq. (4.4.21) in Eq. (4.4.27) 




(4.4.28) 


where the coefficient, B n may vary throughout the vaporization process as 
described in Sec. 4.5. 

(6) Air Properties Required by the Model 


Eight air properties are used in the model. Five of these are constants: air 

molecular weight (ty a ) , critical pressure (P Cfl ) , critical temperature (T c ), and the 
force constants [o a , (e/k) a ], The other three parameters are molecular viscosity 
(u a ), thermal conductivity (k a ) and heat capacity (c Pa ) which are determined from 
polynomial expressions having the form of Eq. (4.4.22). It may be noted that a greater 
level of sophistication in accounting for the temperature dependence of the physical 
properties of the fluid is used in the PTRAK code than in the ADD code. The reasons 
for chis are that the mean temperature (T m ) in the droplet film may be substantially 
different from the air temperature (T a ) and that the mean temperature may vary— due 
to changes in the liquid temperature— even if the air temperature is uniform 
throughout the flow field. Both the air molecular viscosity and thermal conductivity 
Increase by nearly a factor of two as the temperature Increases from 400K to 900K. 
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4.5 Evaluation of Dlatillata Fuel Ef facta 

The vaporization model developed as part of the PTRAK computer program is 
applicable to both distillate and pure fuels. For the present discussion, a dis- 
tillate fuel is a liquid mixture consisting of two or more compounds that may 
differ chemically and physically. The most important property of the fuel for 
determining the vaporization rate is the vapor pressure because of its influence 
on the fuel vapor concentration gradient, the physical basis for the ic 

model. If it is assumed that at any instant, the droplet-air system is in thermodynamic 
equilibrium, then the vapor pressure of the liquid is equal to the partial pressure 
(and, hence, molar concentration) of the fuel vapor at the surface of the droplet. 

For a pure substance, such as water or a single hydrocarbon compound, the 
vapor pressure is a function only of the temperature of the liquid. Thus, as a 
liquid of this type vaporizes, its vapor pressure will vary only if its temperature 
is varying. For a distillate fuel, the vaporization situation is more complicated 
since the vapor pressure is determined by the contributions of each vaporizing 
constituent. To illustrate this complexity, assume that the droplet temperature can 
somehow remain constant. Then the vapor pressure of the distillate liquid tends 
to decrease with time as the more volatile (higher vapor pressure) components 
vaporize so that the vapor pressure of a distillate fuel is a function of both 
droplet temperature and composition. 

The dependence of the vapor pressure of a distillate fuel on temperature and 
composition at a pressure of one atomsphere is given by the standard distillation 
curve such as the curve for a Jet-A type of fuel given in Fig. 4.2. It is seen in 
Fig. 4.2 that the fuel temperature must increase as the more volatile components 
are vaporized at constant pressure. Similarly, another constant value o pressure 
could be used to obtain a distillation curve which would be qualitatively simila 
to the curve presented in Fig. 4.2. Thus, the distillation curve is simply an 

isobar of a pressure-temperature-composition surface. 

Because both the droplet temperature and composition must be known at any 
instant to determine the vapor pressure of the droplet, it is necessary that eac 
droplet in the system be tracked individually. In the PTRAK code, each class of drop- 
lets (a class has homogeneous properties in this context) is treated separate y. 

Thus, in principle the composition of each droplet is known at every instant from 
the distillation curve and the percent of initial mass evaporated. (For droplet 
classes undergoing coalescence, specification of the droplet history requires 

further assumptions.) 

A distillate fuel contains many constituents. It is not practical from a 
computational standpoint to keep track of each component and calculate ^8 con- 
tribution to the vapor pressure. The method used in the vaporization model is based 
upon the work of Cox (Ref. 4.9) and is a standard procedure used in chemical engineering 
practice (e.g.. Ref. 4.10). In Cox’s method it is assumed that, at any instant 
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during the vaporizing process, the behavior of the distillate fuel is the 
same as that of some pure substance; that is. Its vapor pressure is a function 
of the liquid temperature only. Different pure substances characterise the 
vaporization process at each instant. It is assumed that the vapor pressure of a 
hydrocarbon distillate fuel can be determined from the typical distillation curve 
for the fuel and the vapor pressure vs. temperature curves for the normal paraffin 
series of hydrocarbons (C n H 2n +2^ * Thi ® f «®Hy of vapor pressure curves is shown 
in Fig. 4.3. It should be noted that Cox used a temperature scale that allowed 
the vapor pressure for water to be plotted as a linear function of the temperature 
axis; this vapor pressure curve can be approximated by: 


In Pf,s = a n * ~ t ~ * Xn \ * T* 


(4.5.1) 


However, for purposes of this analysis it is more convenient and sufficiently 
accurate to approximate Eq. (4.5.1) by: 


>" Pf,s s a n + TT 


(4.5.2) 


where the subscript n refers to the carbon number of the paraffin which serves as 
the parameter specifying the composition of the distillate at any instant of time. 

Cox's procedure operates as follows. At a particular instant during the 
vaporizing process, the percent of fuel that has vaporized is calculated for each 
droplet. 


(/>i <>!>,„ 


p « = ,0 ° 

[ ipi 0l^=O 


(4.5.3) 


(This calculation requires that the computer program has stored the initial values 
of droplet diameter and density.) The distillation curve (e.g., Fig. 4.2) is entered 
with the fraction of fuel vaporized as a parameter and the distillation temperature, 
T 0 , is found. This is the temperature that is just sufficient to produce a vapor 
pressure of one atmosphere at the specified fraction of fuel vaporized. Of course, 
the droplet temperature is not necessarily equal to the distillation temperature 
nor is the droplet always vaporizing into surroundings maintained at one atmosphere. 
To find the correct vapor pressure, one locates the paraffin curve in the Cox chart 
that passes through the point at a pressure of one atmosphere and the instantaneous 
value of distillation temperature. It is then assumed that this vapor pressure 
curve characterizes the vaporization process at the specified Instant. The instan- 
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taneous value of vapor pressure Is found by moving along the vapor pressure curve to 
a temperature equal to the current value of droplet temperature. The vapor pressure 
is then used to calculate a new value of vaporization rate, droplet temperature, etc. 
(see Section 4.4) and then a new value of fraction of fuel vaporized. The entire 
process is repeated until all of the fuel has vaporized. 


Implementation of the Cox procedure in the PTRAK code is carried out by using 

the form of the vapor pressure curves given by Eq. (4. 5. 2) and inputting two sets of 

values of pressure and temperature for each curve in the family. From these data, 

a and 6 (Eq. 4.5.2) can be calculated — it is assumed that a and 6 vary 
n n 

continuously with some parameter such as the carbon number. Thus, the Cox 
chart is established by two isobars, one of which is taken to be atmospheric 
! pressure (i.e., the pressure used to determine the distillation curve): 


T l , • F, (n , l otm) 


(4.5.4) 


= F2( n i p2^ 


(4.5.5) 



where is any pressure other than atmospheric pressure. For example, 

Eq. (4.5.4) is solved for n with T^ equal to the instantaneous value of 
distillation temperature. Next, Eq. (4.5.5) is solved for T L2 at this value 
of n. Then, Eq. (4.5.2) is applied twice (once at T^ and p^ equal to 1 atm, 

again at T L2 and p 2 ) so that c* n and 8 n are determined. 

Of course, data for any group of compounds may be used to generate a Cox chart 

if that group is appropriate to the multi-component fuel being analyzed. The only 

restriction on the method of utilizing the Cox procedure that has been incorporated 
into the PTRAK computer program is that the two isobars must be continuous functions 
of a variable that describes the group of compounds in some fashion. Note that the 
Cox chart for a pure substance is simply a single vapor pressure curve in the form 
of Eq. 4.5.2. To avoid the use of additional input options, this is accomplished 
simply by inputting data for Eqs. (4.5.4) and (4.5.5) as described in Volume II of 
his report and supplying a distillation curve of the form shown in Fig. 4.4. 
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4.6 High Pressure Effect* 

The model for fuel droplet vaporization and heating Is designed to be applicable 
to a wide range of gas temperature (400 to 900 K) and gaa pressure (3 to 40 atm) . 

As the pressure of a gas Increases (while its temperature remains constant) , the 
density calculated by the ideal equation of state 

p * P^/(R 0 T) (4.6.1) 

begins to differ from the actual density. Futhermore , as the pressure increases t 
the temperature of a vaporizing fuel droplet increases (due to increased convective 
heat transfer) and may approach the critical temperature of the fuel (Ref. 4.11). 

At the critical point, there is no distinction between the liquid and vapor phases. 
These two effects of increasing ambient pressure may lead to significant differences 
between experimentally determined vaporization rates and rates estimated using simp- 
ler droplet vaporization models such as incorporated in the PTRAK code. TVo approaches 
for modeling the droplet heating and vaporization process at high ambient pressure 
were evaluated for use in this analysis. In the first approach (and the one ulti- 
mately incorporated into the PTRAK computer program), mass transfer models are de- 
rived from the diffusion equation: 



which is solved after making simplifying assumptions such as steady flow, spherical 
symmetry, and ideal gas behavior (e.g.. Ref. 4.2 to 4.4). An analogous energy equa- 
tion is also used to derive heat transfer models. The effect of forced convection 
is introduced by modifying the appropriate transfer coefficient by application of 
experimental data. Pressure effects are included either by using thermal and trans- 
port coefficients that may have a pressure dependence or by obtaining solutions to 
Eq. (4.6.2) in terms of partial pressures or concentrations; however, the latter 

requires that the gas-phase behave ideally. 

A second and more rigorous approach begins with the use of the set of differen- 
tial equations applicable to both the liquid and gas (air plus fuel vapor) phase 
(e.g., Ref. 4.12 and 4.13). This set of equations includes an equation of state 
suitable for use at high pressure. It was not within the scope of the present effort 
to modify the computer program used to calculate the air flow field (the ADD code). 
However, since the ADD code uses the ideal equation of state (Eq. 4.6.1), it 
was necessary to assess the departure of air from ideal behavior at elevated pres- 
sures. To do this, the density of air was calculated by both Eq. (4.6.1) and the 
Redlich-Kwong equation (Ref. 4.14). 


^ T ,/2 (v*b)v 


(v-b) 


R 0 T 

~yT 


(4.6.3) 
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where v - 1/ p and a and b are constants depending on only the critical properties 
of the gas (at least for simple molecules). Use of the Redlich-Kwong equation .is 
justified by Shah and Rtodos (Ref. 4.15) who compared various equations of state 
for both argon and n-heptane with experimental data and determined that the Redlich- 
Kwong equation gave superior results. The discrepancies were largest for n-heptane. 
Indeed, it appears that as the molecular structure becomes more complicated, cor- 
rections must be applied to the "constants" a and b (Ref. 4.16). In the present 
study, Eq . (4.6.1) and (4.6.3) were used to calculate the density of the air over the 
range of temperature and pressure cited above and the results agreed within 1.2 per- 
cent of each other. This good agreement is not surprising since the air temperatures 
used (400 to 900 K) are much greater than the critical temperature of air (about 
133 K). At temperatures greatly in excess of the critical temperature, the molecules 
are energetic enough that all gases behave as dilute (and, hence, ideal) gases. 
Furthermore, the Redlich— Kwong equation was developed for gases above the critical 
temperature. 

The good agreement between the ideal and Redlich-Kwong equations of state means 
that it is not necessary to modify the ADD code to incorporate the effects of high 
pressure. However, the gas in the film surrounding the droplet is not air but is 
a mixture of air and fuel vapor (and the mean film temperature may be below the fuel 
critical temperature) . It is conditions in this film which control the vaporization 
rate. If the film consists entirely of a hydrocarbon fuel vapor (such as JP-4 whose 
critical temperature is about 590 K and whose critical pressure is about 31 atm) 
then significant discrepancies between the densities calculated using Eq. (4.6.1) and 
Eq. (4.6.3) result. At 400 K the density calculated from Eq. (4.6.3) is about ten times 
greater than the density calculated from Eq. (4.6.1) for JP-4 vapor. (Even the 
applicability of the Redlich-Kwong equation at such low subcritical temperatures is 
questionable based upon the work reported in Ref. 4.15.) 

The only position in the vapor film at which the vapor could constitute 100 

percent of the mass is at the droplet surface and this can only occur if the droplet 

is at the boiling point of the fuel at the local pressure. Experience at UTRC in 
applying vaporization analyses to problems of fuel-air mixing in gas turbine engines 
indicates that a droplet is not likely to reach its boiling point. A droplet reaches 
a steady-state temperature determined by the maximum vaporization rate compatible 
with the imposed heat flux (the wet-bulb temperature), but the droplet vapor pressure 
at this temperature is considerably less than the compressor discharge pressure which, 

for modern engines, is in the neighborhood of the critical pressure of fuels of in- 

terest. It should be noted that there are conditions where extremely high fuel vapor 
pressure can exist. Calculations by Wieber (Ref. 4.11) indicate that in spray com- 
bustion the droplet temperature can approach the critical temperature (which is 
obviously always at least as great as the boiling point); however, these results are 
based on rocket combustion chamber conditions (temperatures about 3000 K and pressures 
to 150 atm) which are considerably in excess of the flow conditions to which the pre- 
sent model is designed to apply. 
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If It Is assumed that the fual vapor concentration In the film aurroundlng the 
droplet la navar more than 50 percent (by miss), then for repreeentatlve condltlona 
the densities calculated using Eq. (4.6.1) and (4.6.3) are always within 5 percent of 
each other. The discrepancy diminishes to less than two percent at high temperatures 
for JP-4 fuel. (Treatment of this mixture of air and fuel vapor using the Redlich- 
Kwong equation was made In accordance with the rules for determining a and b de- 
scribed in Ref. 4.17). 


As a consequence of these calculations, it is apparent that the use of a mod- 
ified equation of state is not required in this application. Faeth and Chanin 
(Ref. 4.18) suggest an approach to the droplet vaporization process applicable to 
intermediate pressures (about one half the critical pressure) which uses the 
ideal equation of state but which includes real, high pressure, thermodynamic 
effects through the use of fugacity coefficients which are used to specify the 
boundary conditions between the liquid and vapor phase of the fuel. The use of 
such a model presumes the availability of fugacity data and other thermal data 
which are currently unavailable for distillate fuels. Rosner and Chang (Ref. 4.19) 
have considered a vaporizing droplet whose temperature is near its critical point. 
The diffusion equation includes the ratio of the liquid to gas phase density. The 
analysis is applicable to a droplet in quiescent air and thus requires a forced 
convection correction. A numerical solution is required. It was judged that 
simplifying the numerical complexity by solving the problem for a droplet in a 
quiescent medium and then applying a semi-empirical correction for forced con- 
vection to avoid solving a coupled mass transfer and boundary layer calculation 
would reduce the accuracy of the overall approach because of the inaccuracies of 
the forced convection correction. Also, the method would consume significant 
machine computation time. For these reasons, this approach was rejected. 


Because none of the alternative approaches offered significant advantages, the 
first approach was selected for use in the development of the droplet vaporization 
and heating models for the PTRAK computer program. Alternative correlations to 
Ranz and Marshall expressions are used for obtaining mass and heat transfer coeffi- 
cients in the high pressure regime. Following the suggestion of Faeth and Chanin 
(Ref. 4.18), these alternative correlations are used for gas pressures that exceed 
one-half the critical pressure of either the fuel or the air. A correlation suggested 
by Canada (Ref. 4.20) for determining the heat trasnfer coefficient has been incor- 
porated into the model: 1/2 1/5 

ho 0.556 Rep Pr 


Nu h = 


= 2 + 


Re 0 Pr 


237 

d.2/3 


(4.6.4) 


It has been assumed that an analogous correlation exits for mass transfer: 


KDT m R 0 


0 556 Re'p 2 Sc /l 


1.237 


(4.6.5) 
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4.7 Droplet Claes end Distribution Functions 

The system of eight equations Eqs. (4.2.1) through (4.2.6) and Eq. (4. 3. 3) and 
(4.3.9) for the eight dependent variables 


r » s 


location in space 


V V*. V s 


velocity in space 
droplet diameter 


(4.7.1) 


T droplet temperature 

uniquely determines the state of a droplet at any instant of time when these eight 
dependent variables are specified as initial conditions. In principle it is possible 
to solve this set of equations for each of the individual droplets. However, the 
number of droplets can be substantially reduced by treating each droplet as though 
it were a cloud of n fuel droplets the mean properties of which are equal to the 
properties of the individual droplet. This section describes the procedures whereby 
the state of the liquid fuel is described by the behavior of a relatively few number 
of individual droplets. 


Phase Space 


Phase space is nine dimensional space with abscissas given by the eight dependent 
variables, r*, V, D, T^ and the ordinate given by the number of droplets n^ in the I c 
class which satisfies the following conditions 


_ AT 
<r t < T 


- AV - AV I 

— <Vl \ 

( (4.7.2) 

-AD n AD I 

— — < D t < — — 1 
2 2 1 

-ATl , Mi I 

— <Tu ~ ' 

At any instant of time, the distribution of nj over phase space uniquely determines 
the state of the cloud of fuel droplets. 

Droplet Class (Lagrangian Sense ) 

A droplet class is defined in the Lagrangian sense by specifying the number of 
droplets n^ which satisfy conditions (4.7.2) at time t ■ 0. This is equivalent to 
specifying the initial conditions for the droplet (4.7.1) and assigning a number n^ 
to each droplet. In the absence of droplet collisions, the number of droplet classes 
in phase space at any given time determines the state of the entire cloud of droplets. 
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It should ba no tad that in this Lagrangian aystam, more than ona droplet claaa nay 
ba at tha sama location In phaaa apaca. Whan droplata collide or shatter, part or 
all of tha droplata of one class may disappear and a corresponding aaount of fuel 
is assigned to another class. Hence tha nuobar density n^ of a given class nay 
change discontinuous ly. 


Distribution Function 

The distribution function distributes the number of droplets over phase space, 
With a small finite number of droplets, it is often convenient to use an integer 
function called the binomial distribution function with probability of 1/2. This 
function is given by: 


f(l,IL)» 


(1L) ! 

I! (IL-I)! 


I XL 

. I s O.IL 


(4.7.3) 


and has the following properties: 


I f(l,IL) = 10 
I«0 


(4.7.4) 


As IL * ®, the binomial function approaches the normal (Gaussian) distribution 
function. As an example, define the diameter axis in phase space by (see Fig. 4.5) 


D x = D 0 ♦ A D • I 


(4.7.5) 


Then the mean diameter and variance of the cloud of droplets is given by 


D = I 0 X = D x 
1*0 


a o s i ynr ad 


(4.7.6) 


(4.7.7) 


An analogous procedure may be used to distribute fuel droplets over the other coor- 
dinates in phase space. The composite distribution function for all of phase space 
Is formed by the product of the binomial distributions for the eight coordinates in 
phase space, l.e.. 


f i * v y w wv f r 


(4.7.8) 


4 


R82-915362-40 


ORIGINAL PAGE 13 
OF POOR QUALITY 


The composite distribution function Fj hss the property that 


IN-I 

I F t = 1.0 

t:0 


(4.7.9) 


where IN-1 is the total number of classes. Fj may be interpreted as the fraction of 
droplets in the I th class. 

Rosin- Rammler Distribution 

The Rosin-Rammler distribution function is one of several distribution 
functions commonly used to characterize the droplet size distribution typically pro- 
duced by gas turbine fuel injectors (e.g.. Ref. 4.21). This distribution function 
may be used as an alternative to the binomial distribution function Eq. (4.7.8). 

The Rosin-Rammler distribution is defined as the fraction of droplets greater than 
size D and is given by 


=e*P 


-W 


(4.7.10) 


where m is an empirical parameter (typically 1.5<m<3.0). Hence from the definition 
of the distribution function (4.6.2) we have 


(4.7.11) 


dR R dD / 0 \ / D \ 

“-do Sdomm TlT/ exp nT) 


and for a finite number of classes, the Rosin-Rammler distribution function is given 
by 


f-U,IL) = 


/ Di\ m f 

\p/ 6Xp [ \ D / J D 

&(t) •’"MoJ T 


(4.7.12) 


Numb er Density 


The number density is defined as the number of droplets per unit time per unit 
volume of phase space; the total number of droplets is related to the fuel flow rate 


. IN-I IN-I 

w, = £ n x m x * n I f x m x 

1*0 1*0 


(4.7.13) 
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Droplet Class (Eulerian Sanaa) 

The equations of notion for the cloud of fuel droplets is solved In the Lagranglan 
sense, i.e., the fuel droplets ere followed in time and the droplet class is defined 
by its initial conditions. Both the air flow equations and the diffusion equation 
are solved in the Eulerian sense; i.e., attention is fined on a point in space and 
the change in flow in time (or point to point in space) is determined. Therefore, 
in order for the fuel droplet equations to interact with the air flow and gas flow, 
droplet classes in the Eulerian sense must be defined. This procedure is necessary 
because more than one class defined in the Lagragian sense may be in the same element 
of volume of phase space* 

A droplet class defined in the Eulerian sense is the number of droplets in a given 
volume of phase space (Eq. 4.7.2) at a given instant of time. Hence the number 
density of droplet classes defined in the Eulerian sense may change with time. 

The solution of the gas diffusion equation requires only two definitions of 
class in the Eulerian sense; the amount of fuel evaporated per unit volume of space 
per unit time, and the amount of fuel striking the wall per unit area per unit time. 

If dV and dA are the differential volume and area respectively, then 


w f C T) = In,m, for < T < 


(4.7.14) 


. v . -0A 

"fw (r)S l 0 !™! for — 


— dA 
r < — 


(4.7.15) 


where the suranation takes place for only those droplets which satisfy the conditions 
stated. In addition, the state of the cloud of fuel droplets may be described by the 
Sauter mean diameter 


WM .I.M , 

SMD = I F t Df / I F t Dx 
1*0 ' 1*0 


(4.7.16) 
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4.8 Droplet Collision Model 

It has been theorised that the rate of vaporization of a fuel spray may be 
affected by droplet coalescence which can occur upon collision oi. two droplets. 
Coalescence would cause a greater nuaber of smaller droplets to be transformed into 
a smaller nunber of larger droplets; because vaporization time is roughly propor-. 
tional to the square of droplet diameter, the spray vaporization rate is diminished 
by coalescence. Whether coalescence occurs in sprays of the type of interest herein 
is in doubt. Experience gained at UTRC by examining holograms of sprays has shown 
no indication of droplet growth as distance from the injector face increases. In 
the case of extremely dense sprays such as those existing near the face of the nozzle 
of a high- thrust, high-pressure ratio gas turbine engine main burner, coalescence 
may be of significance in the spray development region. However, it is doubtful that 
experimental techniques for acquiring data in regions of extremely high droplet num- 
ber density will be available in the near future. Experimental data do indicate that 
coalescence does not occur for all collisions based upon unreported experiments con- 
ducted at UTRC. In these experiments, droplets of alcohol intersected at right 
angles droplets of fuel oil. The incident droplets had the same diameter (about 100 
to 200 microns). After collision, the droplets coalesced into an unstable configura- 
tion which shattered into two large droplets of approximately equal size and — in 
some cases — one or two very small droplets. (It should be noted that the study o 
droplet coalescence was not the purpose of these experiments.) 

If coalescence is to occur, some energy must be dissipated during the colli- 
sion to form a stable, larger droplet. For droplets of water falling through a 
mist of droplets of the same diameter, Swinbank (Ref. 4.22) has estimated that 
the relative kinetic energy of the colliding droplets is less than one percent 
of the required surface energy for droplets with diameters of 250 microns. 

(These energies would be somewhat less for hydrocarbon fuels because the surface 
tension is less for these substances than it is for water.) Swinbank conducted 
experiments with various fluids and observed no coalescence but his droplet 
diameters were less ^lian 4 microns. 

°“rrr is 

droplet size" increases until there is equilibrium between the rate of coalescence and 
the rate of shattering of the subsequent coalesced mass (Ref. 4.23). 

A review of several papers (Ref. 4.24) stated that "no observed case of a droplet 

In ... coalesce . ' . , th r uot<is collisions between droplets of about 

r~7l» rr “'Lu™ collisions between droplet, o, widel, 
differing sizes do result in coalescence. 


I \ 
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It appears that some of the uncertainty concerning whether droplets coalesce 
Is due to problems of definition. Some workers (such as cited in Ref. 4.24) are 
Interested in collision efficiency; that is, the probability that a droplet will 
collide with other droplets in a given volume. Implicit in this approach is the 
notion that each collision results in coalescence. Yet the previous discussion in- 
dicates that not all collisions result in coalescence. Some Russian workers have 
performed experiments in which a large droplet (400 ym) traverses a cloud of smaller 
droplets (12-13 ym) which have been doped with a special dye (Ref. 4.25). If the 
large droplets were spaced far enough apart so as not to interfere with each other, 
the collection efficiency was about 90 percent; if the large droplets were close 
enough to cause disturbances in the cloud that did not dampen sufficiently prior to 
the traverse by the next large droplet, the collection efficiency was about 40 per- 
cent. A collection efficiency of 100 percent corresponds to a concentration of rye 
in the collected fluid indicating the coalescence of one large droplet and one 
small droplet. 

On the other hand, the experiments in Ref. 4.26 assume that collision efficiency 
and coalescence are highly correlated. In these tests, the droplets were charged 
electrically to simulate conditions occurring in rain cloud formation. Collision 
efficiency was highest for droplets of equal size and equal (but opposite) charge. 

For any charge, the collision efficiency was highest for droplets of equal size. 

Thus, electrical charge appears to have an important effect on coalescence. 

The droplet collision and coalescence problem is too poorly understood and far 
too complicated to warrant anything other than the simplest model. The model pre- 
sented herein is developed in the spirit of Swinbank's approach (Ref. 4.22). 


Assumptions 

In the collision model a number of assumptions are made in order to reduce the 
number of combinations of collisions which are considered and to avoid the definitions 
of new classes. These assumptions are: 


1. Only binary collisions are considered because binary collision are the most 
numerous . 

2. Only collisions between two different classes are considered since it is 
assumed that all droplets of a given class move with the same velocity. 

3. Only collisions between the two classes with the largest number density 
are considered because it is assuned that they are the most numerous. 

4. Only collisions between classes in the same element of volume are 
considered. 


5. If a collision occurs a certain fraction of droplets will coalesce, a cer- 
tain fraction will rebound elastically, and the remainder will proceed 
unaffected. 
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6. To prevent the proliferation of classes, no new classes are created. 
Rather, the properties of the two classes involved in the collision 
are revised. 


Number of Collisions 

In a given volume of space AV, let the Jth class be the class with the largest 
number density and the Kth class be the class with the next largest number density. 
Then the distance between droplets in the Jth class is given by 


/a \ ,/3 

•<$- 


(4.8.1) 


and the number of collisions of class J with class K is 


n COllJK 


T7°k 


(4.8.2) 


Likewise for the Kth class 


= Ur)-°‘ 


(4.8.3) 


n co»IKJ = t K n j 


(4.8.4) 


Then for binary collisions 


n coll = (n C oiijK. n collKj) 


(4.8.5) 


Pr obability of Coalescence 

The probability of coalescence is based on the suggestion of Swinbank (Ref. 4.22) 
and is assumed to be proportional to the kinetic energy of the collision. The mag- 
nitude of the relative velocity of the collision between the Jth and Kth classes is 
given by 


-v K 
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If two droplets coalesce, mass Is conserved and the new droplet 


m * mj ♦ m K 


Is given by 

(4.8.7) 


Likewise the resulting droplet tenperature is given by 


T, = 


ndj C^ L jT|_j 


IC''PUt T lK 


m j p PL.J + m K p PLK 


(4.8.8) 


and the density 


Pl s PlIT l ) 


(4.8.9) 


Thus the coalesced droplet will have a diameter 

Plk P lk 


3 1/3 

|PtJ P LJ * Pl»K P LK ^ 


(4.8.10) 


The work required to expand the Jth and Kth droplet to D is equal to the work re- 
quired to overcome the surface tension force and is given by 


Wj * Sir(0 2 -Dj ) 


(4.8.11) 


W K = Sir(0 2 - 0 2 ) 


(4.8.12) 


Hence the total work which must be dissipated to produce the coalesced droplet is 

»c = WH**V , « (4.8.13) 

Since this energy is dissipated by the collision process, then the probability of 
droplet coalescence is a function of the kinetic energy of the colliding droplets 
due to their relative motion. 


* m j V JK 


/* 


(4.8.14) 




(4.8.15) 
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Hence the kinetic energy available ia 


«C* MIN(«j,«k) 


(4.8.16) 


The probability of coalescence is then defined as that fraction of the total 
number of colliding droplets that undergo collisions energetic enough to dissipate 
the required surface energy: 


p c * C «c /W c * 1 


(4.8.17) 


where C is an empirically determined constant. Hence the number of droplets that 
collide and coalesce is given by 

"c s P C n COll (4.8.18) 

and the number of droplets that collide and rebound elastically is given by 

n RB s ( l ~ Pq ) ^coll (4.8.19) 

The nunber of dorplets not involved in a collision is (n - n coll ) from the Jth class 
and (n K - n col i) from the Kth class. The number of droplets that existed in the Jth 
and Kth classes prior to the collision is given by 


n, = nj ♦ n K 


(4.8.20) 


After the collision, the number of droplets is given by 


n 2 “ P c n c oll 
+ ^“ P c^ n coll 
+ (l-P c ) n coll 
+ (nj-ticoil) 

+ ( n lT n coll ) 

“ "j+Wcoll 


(number coalesced) 

(number elastically rebounding — Jth class) 
(number elastically rebounding - Kth class) 


(number not involved- Jth class) 
(number not involved-Kth class) 


(4.8.21) 


Hence the total number of droplets is reduced by the amount 

n 2 = n i p c n coll 


(4.8.22) 
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Elastic Collisions 


SSS35W 


For nonrotating droplets moving In a three-dimensional space (eae Fig. 4.6), 
the dynamics of a collision can be reduced to a one dimensional problem by resolving 
forces along the line of centers. If V JX and V are the velocities before a col- 
lision and Vj- and V r2 are the velocities after a collision, then along the vector 


I1D1UU ttUU V J2 

t (see Fig. 5.6) 


m K v CK , (i-e)+(mj-em K ) v CJI 


mj + m K 


(4.8.23) 


mjVcji d-e) -Monj-em^VcKi 

mj + m K 


(4.8.24) 


where e is the coefficient of restitution (Ref. 4.27). For perfectly elastic col- 
lisions e ■ 1 and for perfectly inelastic collisions e * 0. 

Class Properties After Collision 

From Eqs. (4.8.20) and (4.8.21) it can be seen that the collision of two classes 
affects conditions in five classes. In order to avoid creating new classes after 
every collision, the Jth and Kth class are redefined with new properties which are 
weighted averages. Since the fuel flow rate of liquid droplets is preserved through 
a collision, 


W IJK = n ji m ji * °KI m Ki = Pc "coll w * (l ' p c)ncoii + (n j _n coiO r ''j 


(4.8.25) 


4-0 - P c >n C0 |i m K 4 -(n K -n C 0 ||)m K 


In the present discussion, it is assumed that the Jth class was the more numerous 
class prior to the collision. All of the coalesced droplets are arbitratily assigned 
to this class. Then the mass and number densities in the Jth and Kth classes can be 
determined from Eqs. (4.8.26) through (4.8.30): 


n j2 m j2 * P C ncc.™ ' Pc>"coll Ocoiilmj, 


(4.8.26) 


n K2 m K2 = ** “ P C * n COll * n jl “ n COlP m j! 


(4.8.27) 
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n J2 * n JI 

n K2 = n Kl“ P C n COll 

°2 “ °l “ P C n c<rfl 


The new distribution functions can be defined by 



fj2 * n J2 / n 2 

(4.8.31) 


*K2 * n K2^n 2 

(4.8.32) 

The remaining properties of the new J class are given by 



P C n C0ll m ^C + U"P(cJ n COll m j| V j2 + (n ji -n coll )m ji v ji 
JZ ~ P c n co ,|7¥? + (l- P c ) n C0 || nrtj, ♦ (rij, -n co iOmji 

(4.8.33) 


• 

Pc^coii™ + (1- Pc)o CO ||mj,C PL + T LJ , + (Oj, -n CO ||)mj,Cp L T LJ i 

(4.8.34) 

T LJ2 = 

P c n c0 „ ?nC PL + (I -P c ) n C olinr)j|CpL -M",,, - n c0 H>n»j,Cp L 

Likewise 

for the Kth class p - o l J ) 

r LJ2 ' L ' LJ2 1 

(4.8.35) 


/ \' /5 



/ 6m ja\ 
°J2 ~ \ Puz ) 

(4.3.36) 


_ !•" ? cl n COll m KI V K2 * * n KI " n COll )m KI V KI 

v k2 = 

(i-P c )n C0 || m«i ♦(n R) - n CO ||)m K | 

(4.8.37) 


(4.8.28) 

(4.8.29) 

(4.8.30) \ 
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(4.8.38) 


(4.8.39) 

(4.8.40) 
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4.9 Droplet Stuttering Model 

Under certain circumstances, the aerodynamic and friction forcea acting on a 
droplet will cause enough Instabilities In the liquid so that a larger droplet 
shatters into smaller droplets. Examination of the literature (Refs. 4.28-4.29) 
indicates that large droplets (*i» 10^ ym diameter) subjected to high aerodynamic 
forces (as produced using shock tubes) are most likely to undergo breakup. It is 
not clear whether the relatively small droplets injected into the low Mach number 
fuel preparation section of a gas turbine can be shattered in this manner. Never- 
theless, a droplet breakup model has been included in the computer program. 

The droplet shattering model used in die analysis is based on a model proposed 
by Wolfe and Andersen (Ref. 4.28). In this model it is assumed that the breakup 
of the liquid droplets is controlled by rate processes; hence, the droplet must be 
subjected to stresses for a period of time before breaking into a finer ensemble of 
droplets. Breakup is assumed to occur because of deformation of the initially spheri- 
cal droplet due to aerodynamic forces (bag breakup) and to stripping of liquid from 
the droplet surface due to friction forces (shear breakup) . The time required for 
droplet breakup to occur is given by: 


*B = 


where : 


[a 2 + bp] ,/2 -a 

I6/i l 


A = 


p Q 0 


and 


B = 2 /p L 


p * k . J Pm c o l Au | 2 -k 2 s/0 


(4.9.1) 

(4.9.2) 

(4.9.3) 

(4.9.4) 


The values of the constants kj^ and k 2 are the empirically determined values 
of Wolfe and Andersen 


Breakup Mode 

k l 

k 2 

Aerodynamic Drag 

0.333 

2.0 

Friction Drag 

0.667 

4.0 


The values of the breakup time are computed at each step in the trajectory 
for each droplet class; the minimum of the aerodynamic or friction breakup time 
is compared with the residence time accumulated by the droplet class subsequent 
to its being subjected to breakup forces. The zero time reference point for 
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breakup time la taken to be that previous Instant of tins at which the calculated 
breakup tine undergoes a large step decrease In magnitude; this would correspond 
to the instant of time at which the droplet was Injected Into the airstream. 


The diameter of the droplets produced after the shattering process (D2) la 
given by an equation suggested by Wolfe and Andersen: 


f *3 Ml s 
L Po p'l 


|AU | 4 


(4.9.5) 


A value of k-j equal to 136 as suggested by Wolfe and Andersen is used in the model. 
The number of droplets in the class undergoing shattering is adjusted by: 


= "(¥) 


(4.9.6) 


so that the mass in this class is conserved. 
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4.10 Droplet-Solid Boundary Interaction Model 


The trajectory calculation used In the PTRAK computer program Is used to 
determine whether a droplet has struck a wall of the fuel preparation section 
passage. A fraction of the incident particles will rebound elastically from the 
wall; the other particles will adhere to the wall and may undergo evaporation. 

The probability of rebound PRBJ for the Jth class is the fraction of incident 
particles that rebound elastically. Two mutually-exclusive probability functions 
are available in the model: 


0.0 < PRBj = C, < 1.0 


(4.10.1) 


0.0 < PRB. = I -C 2 (1/ W ). <1.0 


(4.10.2) 


where and C 2 are input constants. The parameter, I/W, is the ratio of droplet 
kinetic energy to deformation energy for the Jth class. 

In class J, the magnitude of the droplet velocity is VJ so that the droplet 
kinetic energy that must be dissipated at the wall is: 


t - 1 ,.2 w °J 

* " 2 ^ LJ 6 


(4.10.3) 


If the droplet (upon striking the wall) flattens into a hemisphere of diameter DJ , 
then conservation of mass requires 


Dj = 2 I/3 0j 


The total surface area of the flattened droplet is 


(4.10.4) 


a' = \ imli ! * 5- = -f -FTtoJ)* 


(4.10.5) 


whereas the droplet surface area prior to deformation is 


Aj = fr Dj 


(4.10.6) 


% 
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so that the work (energy) of deformation la 


Thus— 


w = S(a'-A) = Siro,? 


ti/wjj «0.4 /DLj q,v*/s 


(4.10.7) 


(4.10.8) 


The I/W ratio may be interpreted to mean that a certain amount of kinetic energy 
must be dissipated on impact in order that sufficient deformation of the droplet 
occur to cause it to stick to the wall. 

For the fraction PRBJ of the droplets that rebound elastically, the droplet 
properties are identical to those of the incident droplets except that the normal 
velocity component is -VJn. The fraction 1-PRBJ of droplets striking the wiTT^main 
thereon and may evaporate. It must be remembered that the overall problem for the 
uel mixer is a steady-state problem. The droplet transients are only important for 
relating droplet lifetime to spatial locations within the mixer flowfield. Thus, 
it is not necessary to consider the rate of fuel evaporation from the wall (as 
described, say, in Ref. 4.31). Rather, it is only necessary to determine how much 
fuel has evaporated at location X on the wall; this amount will serve as a boundary 
condition in the fuel vapor diffusion equation. 

Any fuel that does not evaporate but remains on the wall as a film will (in 
time) reach a steady-state temperature equal to the wall temperature, Tw(x). 

To determine the fate of this fuel, the percentage of fuel evaporating from the 
wall is guessed and the distillation temperature is found from the distillation 
curve. The Cox chart is then used to calculate the vapor pressure at Tw(x) 

(see Sec. 4.5). If the vapor pressure is higher than the local ambient pressure, 
the guessed value of fuel evaporated is too high; if the vapor pressure is less 
than the local pressure, the guessed value is too low. The guessed value is 
adjusted until the vapor pressure and local pressure agree. 
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4.11 Numerical Methods 

The set of equations Eqs. (4.2.1) through 4.2.6) and Eqs. (4.3.3) and 
(4.3.9) are solved using a predictor corrector method which is second order 
accurate in At. This algorithm may be described in functional form as follows: 


^♦i 3 + ^ I r m » » \m » ®cn I ^t 


(4.11.1) 


_ _ »4| 

r l4i 3 r j ¥ 


(4.11.2) 


\i*i ~ \i + ® 


(4.11.3) 


V + I — • — if x) 

~ + R l r m*^m »\«n*0in) At 


(4.11.4) 


3 7 <*T-V * 


(4.11.5) 


1 -*► ♦! — . 

2 ^ r l ♦* + r l * 


(4.11.6) 


I *♦! 

* 7 < \ t - 1 LI> 


(4.11.7) 


*+1 I * +1 

Dm 3 I0j4| + I 

AS/V 


(4.11.8) 


(4.11.9) 


The initial guess for the variables at station 1+1 are the corresponding values 
at station I. Then the differential equations are evaluated at the mid point. 
The iteration continues until 


At ~ At 


(4.11.10) 
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4.13 List of Symbols 

Droplet cross-sectional area (cm^> 

Droplet surface area (cm^) 

Mass fraction of gaseous fuel (dimensionless) 

Drag coefficient of droplet (dimensionless) 

Specific heat at constant pressure (cal/gm/°K) 
Specific heat at consent volume (cal/gm/°K) 

Droplet diameter (cm) 

Mass diffusivity (cm^/sec) 

Coefficient of restitution (dimensionless) 

Droplet cloud distribution function (dimensionless) 
Droplet heat transfer coefficient (cal/cm^/°K/sec) 
Droplet kinetic energy (gm enr/sec*') 

Thermal conductivity (cal/cm/°K/sec) 

Mass transfer coefficient (gm/cm2/atm/ 0 K) 

Distance between droplets (cm) 

Droplet mass (gm) 

Molecular weight (dimensionless) 

Norma 1 coord inate (d imensionless) 

Droplet total number density (l/sec) 

Droplet number density Ith class (1/scc) 

Number density of collisions (l/sec) 

Number density of elastic collisions (l/sec) 
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List of Symbols (Coat'd) 

Nusselt number for heat transfer (dimensionless) 

Nusselt number for mass transfer (dimensionless) 

Probability of coalescence (dimensionless) 

Or, critical pressure (atm) 

Percent of liquid evaporated (dimensionless) 

Pressure (atm) 

Vapor pressure (atm) 

Prandtl number c p u/k (dimensionless) 

Probability of rebound (dimensionless) 

Heat transfer rate (cal/sec) 

Radius from axis of symmetry (cm) 

Position vector for droplet (cm) 

Droplet Reynolds number (dimensionless) 

Universal gas constant (82.0575 cm^.atm/mole/°K 

Rosin-Rammler function (dimensionless) 

Streamwise coordinate (dimensionless) 

Surface tension (dyne/cm) 

Schmidt number (u/o/4) (dimensionless) 

Sauter mean diameter (cm) 

Time (sec) 

Droplet breakup time (sec) 

Temperature (°K) 

Critical temperature (°K) 
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List of Symbols (Coat'd) 

Distillation temperature (°K) 

Air velocity components (cm/sec) 

Volume (cm^) 

Droplet velocity components (cm/sec) 
Metric scale coefficient (1/cm) 

Droplet velocity vector (cm/sec) 
Vaporization rate gm/sec 
Evaporation rate per volume (gm/cm^/sec) 
Evaporation rate per area (gm/cm^/sec) 
Work done by surface tension (dyn cm) 
Liquid fuel flow rate (gm/sec) 

Mole fraction (dimensionless) 

Blowing parameter (dimensionless) 

Mass transfer parameter (dimensionless) 
Cox chart parameters for vapor pressure 
Force constant (°K) 

Heat of vaporization (cal/qm) 

Molecular viscosity (gm/cm/sec) 

Density (gm/cm^) 

Force constant (X) 

Variance in binomial distribution (pm) 
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Subscripts 


f 

L 


s 


0 

1 

2 

I,J,K 



List of Symbols (Cont'd) 
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Air properties 

Gaseous fuel properties 

Liquid fuel properties 

Mean (air + fuel vapor film) properties 

Conditions at droplet surface 

Initial conditions 
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* r? 


Conditions before collision 


Conditions after collision 


I, J, K th class 


Conditions far from the droplet surface (at outer edge of film) 
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TABLE 4.1 

COLLISION INTEGRALS FOR CALCULATING BINARY DIFFUSION CO-EFFICIENTS 


T* 

nm 

T* 

8 < T *> 

0 30 

2.662 

2 7 

0.9770 

0 35 

2 476 

2.8 

09672 

0 40 

2.318 

2 9 

09576 

045 

2 184 

30 

0.9490 

050 

2.066 

3 1 

09406 

0.55 

1 966 

32 

C 9328 

060 

1 877 

3.3 

09256 

065 

1 798 

3 4 

0 9186 

0.70 

1 729 

3.5 

0.9120 

0.75 

1 667 

3.6 

09058 

080 

1 612 

3,7 

08998 

085 

1 562 

3.8 

08942 

0.90 

1517 

3.9 

08888 

095 

1.476 

4 0 

08836 

1 00 

1 439 

4.1 

08788 

1 05 

1 406 

4.2 

08740 

1 10 

1 375 

4 3 

0 8694 

1.15 

1 346 

4.4 

08652 

1 20 

1 320 

4.5 

0 8610 

1 25 

1 296 

4 6 

08568 

1 30 

1 273 

4 7 

08530 

1 35 

1 253 

4 8 

0 8492 

1 40 

1.233 

4 9 

0.8456 

1 45 

1 215 

5 

08422 

1 50 

1.198 

6 

08124 

1 55 

1 182 

7 

07896 

1 60 

1 1 67 

8 

07712 

1 65 

1 153 

9 

0 7556 

1 70 

1 140 

10 

0 7424 

1 75 

1 128 

20 

06640 

1 80 

1116 

30 

0 6232 

1 85 

1 105 

40 

0 5960 

1 90 

1 094 

50 

0 5756 

1 95 

1 084 

60 

05596 

2 00 

1 075 

70 

0 5464 

2 10 

1 057 

80 

0 5352 

2 20 

1 041 

90 

0 5256 

2 30 

1 026 

100 

0 5170 

2 40 

1 012 

200 

0 4644 

2 50 

0 9996 

300 

0 4360 

2 60 

0 9878 

400 

0 4170 
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FIG. 4.1 


GEOMETRY OF DROPLET VAPORIZATION MODEL 
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5.0 ANALYSIS OP DIFFUSION OF FUEL VAPOR 


5 . 1 General Approach 

Since the mass fraction of fuel vapor is snail (lean equivalence ratio), the 
diffusion equations reduce to one equation for the mass fraction of fuel vapor 
diffusing Into a known gas (air) flow field. The thermodynamic properties of 
the fuel vapor and air mixture are taken to he the thermodynamic properties of 
air. It is further assumed that the cloud of liquid fuel droplets does not 
interact with the fuel vapor other than to produce a source term In the diffusion 
equation due to evaporation and a source of diffusion flux at the wall due to 
evaporation of fuel from the wall. Since the flow field Is turbulent, the 
turbulent diffusion coefficient Is assumed to be proportional to the effective 
turbulent viscosity which is calculated In the solution of the air flow field. 
Under these conditions the diffusion equat ton reduces to a Itnear second order 
parahol ic partial differential equation which can be solved bv a forward marching 
integration technique. The Vapor Diffusion IVAPDIF) code Is used to solve 
this diffusion equation. 

vmho 1 s usod In this soot ion ,;ro dot i nod on pa«e 7 \. 
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5.2 Fuel Vapor Diffusion Equation 

If diffusion in the streamwise flow direction is neglected, then the diffusion 
equation in orthogonal streamline coordinates is given by: 

* s dS P dn r d+ r )dn[s c dn | d4|Scrv*df)j f 


J n*“ 


bt 

Sc 


ac 

an 


W 


fw 


(5.2.2) 


Since the gas flow field (p , U s , U n , U^) is known, from the ADD code cal- 
culation and w ^ and w^are known from PTRAK calculation, Eq. (5.2.1) is a linear 
second order partial differential equation in the mass fraction of fuel C which 
can be solved by forward marching techniques with the boundary condition J n 
at the wall given, and the initial conditions C (n, $) prescribed. 
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5.3 Numerical Methods 
The computational mesh Is defined by 


s z * As(l-i ) 

(5.3.1) 

n j * n(J) 

(5.3.2) 

♦k s^IK-I) 

(5.3.3) 

An, * n(j+i) - n(j) 

(5.3.4) 


where it is noted that the computational mesh in the n direction is not uniformly 
spaced. For a nonuniform mesh the difference formulae are given by 


* [W“ (, ‘ r2) V’ r2 W]/ d « (5.3.5) 

-~5 = jtiai.K ”( , * r ) fi,K vr fj-l,it] /&2 (5.3.6) 

r * Arij Mnj.| (5.3.7) 

d|* An, ♦r* An 4 _, (5.3.8) 

d 2 * (&nj + r AnJ., >/ 2 (5.3.9) 

If * ( C - 


where f is any variable. These difference formulae are also applicable to a uniform 
grid where r - 1 in the direction. 


Equation (5.2.1) is solved using the method of point successive over relax- 
ation (see Roache (Ref. 5.1)). Let the difference formula for mass fraction 

C be given by 



hL 
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dC | 

dn * 1 

9i.it ~ r*<s_r,«]/ d ( 

(5.3.11) 


d 2 C 

77 = 

f ¥ ¥ + \ 1 / 

| c jfi,K — Cl -r ) Qi.k**" r Cj-i.kJ/ da 

(5.3.12) 


Then the difference formulae Eq. (5.3.5) through Eq. (5.3.12) can be applied to 
Eq. (5.2.1) and solved for as the unknown in terms of the value of C at 

neighboring points. The v+1 guess for is then 

C J,K = Cj,k “Cjk) (5.3.13) 


The relaxation factor R as determined by Garabedian (Ref. 5.2) is given by 


<R = 2/(1 + 3.014 h/A l/2 ) (5.3.14) 

where 

A = njL^Kl (5.3.16) 


Equation (5.3.13) is applied to all Interior points. At the boundary, the one- 
sided difference formula is applied to Eq. (5.2.2). As an example, for J - 1 


* [~ C 3,K -f ( r ♦ 2)C^']y , d, 

\J+1 

which can be solved for C 

J,K 

Iteration continues until 



(5.3.17) 


♦! 


- C* 


C v 


< t 


(5.3.18) 
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5.5 List of Symbols 

Mass fraction of fuel vapor (dimensionless) 

Arbitrary variable 

Normal coordinate (dimensionless) 

Radius to axis of symmetry (cm) 

Streamwise coordinate (dimensionless) 

Schmidt number (dimensionless) 

Air velocity components 

Evaporation rate per unit volume (gm/cmVsec) 
Evaporation rate per unit area (gm/cn 2 /sec) 
Metric scale coefficient (1/cm) 

Tangential coordinate (dimensionless) 
Effective turbulent viscosity (gm/cm/sec) 
Density (gm/cm 3 ) 
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6.0 COMPUTER MODEL CALIBRATION 


6.1 Calibration of Fuel Droplet Model 

in Ir^T W6re t0 6 acqulred at conditions simulating those occurring’ 

Lfrrrnr 8 ^, ts and ue " to *■ -«*«*«* 

measures of the initial conditions required by the calculating procedures. 

air f^irr mentS W6re COnducted ln an 8 -89 x 8.89 cm square channel in which 
d ^ at p Pr ? SSUreS ^ tem P erat ures representative of compressor discharge 
conditions Fuel was injected through an orifice located on tL d^st^e^ side 
of a circular tube at its mid-span point. The tube traversed the ™tlTf It 

in C L an vi ci^ro^hehoriflce? SymnetrlCal 3ir fo11 ov « lts -tire length except 

to UTRC Jet NASA y nhf- a J rC r ft ^ tUrbinC fuel was used in the experiments reported 
1™°'. obtalned some Properties of the batch of Jet-A fuel used in these 

to ^Rc l r k ln TabU h* 1 * Fuel P ro P ert ies not supplied by NASA 

mated by UTR^from Perf0nned *“ Ta8k ” “ d ™ k 111 We ~ e8ti ‘ 

„ m f Dat l W T : eCOrded at three 8x131 positions that were located 7.5 15 and 30 

Of *7 ' l T C ° r • Sta8Mtl ° n P rM8 “" “ d 

o air were measured at several points along bisectors of the duct both paral- 
lel and perpendicular to the axis of the fuel injector strut. Laser £opp£r 
ve ocimetry (LDV) was used to determine the axial and radial components ^f air 

si Z e C d^str^!t^ al t and COmp ° nentS of dr °P let velocity, and the fuel droplet 

ze distribution at up to sixteen locations in the measuring half-plane at each 

A spulover t€chni,ue - "* ed *• *»« fo-i i:i„: r ‘« lon 

*,d.l T ?!nK ta ,r tUal ^, aC ' ,Ulred ln tha e *P erl ~" c « proved to be Inadequate for 
odel calibration. Only one set of data «as provided to UTRC, and because of 

arge droplet sizes involved and hence very low levels of vaporization the 
data were of limited use in calibrating the model. Data were acJuireJat a pres- 
sure of 5 atm and at a temperature of 650K. The axial velocity of the air was 

D^obl^r 1 til Jh m A SeC ' ThC radial veloclt y of the alr not reported because of 
Ih t a U h , ^ data reductlon equipment. Data for the fuel droplets indicated 
that the droplet axial velocity was typically 30 m/sec at the 7.5 cm measuring 

station. Duplet radial velocities were typically a few meters per second at 
this point. The spray velocity and mean particle size distributions were not 


74 


R82-915362-40 


symmetric about the mid-plane of the duct parallel to the fuel-injector strut 

perhaps due to a misalignment of the fuel injector strut. The data indicated i j 

that the mean droplet size 278 urn at the 7.5 cm measuring station) was essen- j. .'jj 

tially unchanged between the 7.5 to 30 cm locations. Fuel vapor concentration 

measurements were not reported to UTRC. j i 

These data could not be used to calibrate the model since essentially no i J 

vaporization occurred between the 7.5 and 30 cm measuring station. The computer 1 

codes were used to verify this result. The geometry of the test section together 
with the air flow conditions were input into the Annular Diffuser Deck (ADD) 

code. The effects of the str it were ignored since it was believed that the strut .j 

would have no measurable impact on the calculated vaporization rates. The drop- ^ 

let axial and radial velocity components together with a mean droplet size and j 

variance were input into the particle tracking (PTRAK) code at each of the twelve j 

locations reported for the 7.5 cm measuring station. The program calculated that : 

less than 3 percent of the fuel vaporized between the 7.5 cm and 30 cm axial \ 

positions. • 

The heat transfer, mass transfer and drag coefficients used for these cal- 
culations were those discussed in Section 4.4 and have been used by many other ] 

workers in this field. The procedure used to determine the vapor pressure of the * , 

fuel (Section 4.5) has not — as far as it is known — been used by other workers pre- 
dicting fuel droplet vaporization behavior; however, the technique is a standard \ 

chemical engineering procedure. While it is comforting that the calculated and 
experimental results agree, it must be recognized that the vaporization and tra- 
jectory calculations can be expected to require calibration at least for flows : 

at moderate pressures. Unfortunately, the test program had to be terminated 
before data could be acquired for calibration of the models used to determine the 
effects of droplet shattering, droplet coalescence, droplet-solid boundary inter- j 

actions (other than simple elastic rebounds) and high gas stream pressure. a 
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6.2 Model Sensitivity Study 

The intent of the model sensitivity study was to assess the behavior of the 
model (the system of three computer programs) when flow conditions, empirical 
factors, and numerical factors were altered about selected mean values. Empirical 
factors, such as the heat and mass transfer coefficients used in the vaporization 
analysis, were to be determined (or verified) as part of the task of model calibra- 
tion. Numerical factors are those parameters that control the mathematical pro- 
gress of the codes such as those used to determine the step-size in forward- 
marching integration procedures. 


There are three computer program used in the model. The gas flow field is 
determined by the Annular Diffuser Deck (ADD) code. The behavior of the fuel 
spray is calculated using the particle tracking (PTRAK) program. The resulting 
distribution of fuel vapor is estimated using the Vapor Diffusion (VAPDIF) deck. 
The results of VAPDIF are determined once the fuel vapor source terms and the 
computational grids established by the ADD and PTRAK codes are specified. The 
ADD code is a well-established code in use at several government and industrial 
sites. The sensitivity study was confined to evaluating the influence of changing 
air or fuel spray conditions on predictions made using the PTRAK code. 

The baseline conditions used in the sensitivity study were essentially those 
of the calibration experiments; air stagnation pressure of 5 atm, air stagnation 
temperature of 6S0K, fuel injection temperature of 310K, and fuel droplet Initial 
velocities of (approximately) 30 m/sec. The droplet velocity varied from point 
to point in the initial calculating plane similar to the experimentally-determined 
variation. The calculations were performed over a distance of 22.5 cm corres- 
ponding to the distance between the 7.5 cm and 30 cm measuring stations used in 
the calibration experiments. 

Effect of Initial Droplet Diameter 


The reported droplet Sauter mean diameter (SMD) for the calibration experiments 
(at the 7.5 cm measuring station) was 278 ym. The model predicted that less than 
three percent of the fuel evaporated over a distance of 22.5 cm. Calculations 
were also carried out for initial droplet Sauter mean diameters of 100 and 50 ym; 
the model predicted that 13 percent and 39 percent, respectively, of the fuel 
would vaporize. All subsequent calculations were performed with an initial SMD 
equal to 50 ym. 


E f f ect of Step-Size 


Increasing the step-size in the particle tracking code by a factor of 2.5 
above the value used in the baseline case had no effect on the calculated amount 
of fuel vaporized (39 percent). Step-size effects are expected to be significant 
only for cases in which the droplets are small and have velocities significantly 
different from the local air velocity. 


76 


R82-915362-40 


Effect of Initial Fuel Temperature 

If the fuel is preheated to 360K, the analysis indicates that 44 percent of 
the fuel vaporizes. This relatively small effect of fuel initial temperature is 
due to the droplets quickly reaching a temperature at which the rate of vaporization 
is just supported by the imposed heat flux. Thereafter, the droplet temperature 
in either case increases slowly. For droplets of pure substances evaporating in 
a uniform gas environment, the droplet temperature may reach a constant (wet-bulb) 
vaLue. For droplets of a distillate fuel (the case analyzed here) in the same gas 
environment, the droplet temperature will never reach a constant value less than 
the ambient temperature, the droplet temperature will increase continuously as 
the more volatile components are evaporated. The calculated fuel droplet histories 
are shown in Fig. 6.1. For cases in which the initial droplet temperature 
transient is substantially longer, the results would differ more significantly. 

Effect of Fuel Volatility 

A more volatile fuel can be expected to vaporize more rapidly than a less 
volatile fuel. Fuel volatility was stimulated by reducing uniformly the fuel 
distillation curve by 30K. The analysis predicted that 35 percent of the fuel 
would vaporize. 

Effect of Ambient Temperature 

The stagnation temperature of the gas was increased from 650K to 750K. The 
air velocity was maintained at the original value. The ADD code was used to 
generate the new flow field. The PTRAK analysis calculated that 63 percent of fuel 
vaporized. This result is reasonable since the heat transfer rate to the droplet 
increases with increasing ambient temperature. 

Effect of Initial Droplet Axial Velocity 

The initial axial velocity of the droplets was reduced to 10 m/sec from 30 
m/sec; the PTRAK code calculated that 48 percent of the fuel would vaporize. Some 
of this increase in fuel vaporization is due to the increased relative velocity 
and, hence Reynolds number, based on droplet relative velocity; the droplet heat 
and mass transfer coefficients used in the analysis (see Section 4.4) are strong 
functions of Reynolds number. Droplet residence times increase when the droplets 
have lower initial velocities, (Fig. 6.2), and the increase in residence time 
augments the effects of increased Reynolds number. For smaller droplets the 
effect of initial droplet velocity would be less important because these particles 
would accelerate more rapidly to the local gas velocity. 

Effect of Ambient Pressure 


As the ambient pressure of the gas Is increased, the droplet Reynolds 
number increases and the heat and mass transfer coefficients increase in magnitude. 
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Thus, an Increase In ambient pressure can be expected to increase the droplet 
vaporization rate* On the other hand, the dreg force on the droplet eleo increases 
because the decrease in drag coefficient is aore then offset by the increased 
dynamic pressure acting on the droplet. Por a droplet with an initial velocity 
less than that of the air, the droplet accelerates sore rapidly as the ambient 
pressure increases and the droplet residence time in the premixing passage de- 
creases. The effect of Increasing ambient pressure on the extent of fuel 


vaporized is the integrated result of an increased vaporisation rate and a de- 
creased residence time. No statement can be made a priori regarding the effect 
of increasing ambient pressure. 


In the case analyzed the air stagnation pressure was Increased from 5 to 
10 atm. The extent of vaporization was reduced slightly from 39 to 38 percent. 
Evidently, the decrease in droplet residence time had a slightly greater influence 
than the increase in vaporization rate. 

Simultaneous Effects of Ambient Pressure and Initial Droplet Velocity 

In this case, the ambient pressure was again increased from 3 to 10 atm 
but the initial droplet axial velocity was decreased from 30 to 10 m/sec. The 
lower initial droplet velocity Increased the residence time in the premixing 
passage and the PTRAK code calculated that 43 percent of the fuel vaporized. 
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6.3 Calibration of Vapor Diffusion Model 

Calibration of the vapor diffusion model under Task II using the VAPDIF code 
was undertaken using experimental data provided to UTRC by NASA-Lewis Research 
Center. These data were obtained by NASA under a separate contract with Solar 
Turbines International. 



The experiments were conducted in the test facility described in Section 6.1 
with methane injected into the flow through a small orifice at the initial 
(Z = 0 cm) station. Methane concentrations were measured over the cross-section 
of the channel at the grid points shown in Fig. 6.3 at three downstream locations 
(Z = 7.5, 15, 30 cm, respectively). The average mass fraction and fuel flow rate 
at each station were determined by integrating the concentration data over the 
cross-sectional grid. These results are presented in Table 6.2 where it will be 
seen that the flow rate of methane increases by 60 percent from the first to 


third measuring station. 


The VAPDIF code was run using the flow conditions presented in Table 6.2. 

For calculation purposes, the flow field was divided into a grid with 25 x 50 
points at each of 38 axial locations. Convergence occurred in about 20 iterations 
with the residual less than 10"^. Computing time was about 20 sec/station. 

Initial methane concentrations were obtained by curve fitting a function to the 
data recorded at the Z = 7.5 cm station. 


The calculated methane concentration (mass fraction) distribution at the 
Z = 30 cm plane is shown in Fig. 6.4. This distribution is similar to the typical 
bell-shaped curve obtained from a point source. This result is not surprising 
since in the experimental program the methane was injected into the flow field 
through a small orifice located on the rig centerline. Comparisons of the cal- 
culated and measured concentration distributions are shown in Fig. 6.5. The 
data along the X = 3.61 cm grid line (see Fig. 6.4) was selected for comparison 
purposes . The results shown in Fig. 6.5 were obtained using a Schmidt number equal 
to 0.10. Only general agreement could be obtained. 


Higher values of Schmidt number produced poorer agreement with the data. 
Generally, the turbulent Schmidt number is of the same order as the molecular 
Schmidt number (Ref. 6.2) which is about 0.8 for methane and air mixtures at 
"he nr “sure and temperature presented in Table 6.2. Since the data indicate a 
lrck of mass conservation and since the Schmidt number varies inversely with the 
, n , ss diffusion coefficient, the VAPDIF computer program could not predict accurately 
the downstream concentration profiles given the concentration profile at the 
7 = 7.5 cn station. 




\ 
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6.5 Tables 
TABLE 6.1 

Properties of Jet-A Fuel Used in Calibration Experiments 


Kinematic Viscosity 
Specific Gravity 
Pour Point 
Flash Point 

Carbon Residue on 10 percent Residum 
Distillation Curve 

Initial Boiling Point 

10 percent distillation 

20 percent distillation 

50 percent distillation 

90 percent distillation 

End point (98.5 percent recovery) 


1.55 cs <3 311K 
0.8096 @ 289K 
211K 
333K 

0.28 percent 


459K 

472K 

479K 

494K 

524K 

560K 
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TABU 6.2 
EXPERIMENTAL DATA 
Nathan* Gaa Diffusion Data 

Air Velocity V - 39.947 n/aac 
Air Density p - 2.7929 Kg/« 3 

Duct Area A - 79.02 M 2 

Weight Flow W a - .8816 Kg/sec 
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FIG. 6.2 


CALCULATED DROPLET RESIDENCE TIMES 


■ DROPLET INITIAL VELOCITY = 30 m/sec 
— — DROPLET INITIAL VELOCITY = 10 m/seo 
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FIG. 6.3 
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COMPARISON OF CALCULATED AND MEASURED METHANE CONCENTRATION 
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7.0 PARAMETRIC STUDIES 


The model , consisting of three compuater programs, was applied to the 
analysis of two configurations for premixing passages operating at conditions 
corresponding to Energy Efficient Engine (E^) cruise conditions. In the analysis 
of both designs, it was assumed that the inlet air stagnation temperature was 
745K, the inlet air stagnation pressure was 11.1 atm, and the inlet air mean axial 
velocity was 54.9 m/sec. The equivalence ratio used in the analysis was different 
for each case. 

Symbols used in this section are defined on page 92. 

7.1 Swirl Tube Premixing Passage 

In the swirl tube design, fuel is injected via a pressure atomizer along the 
centerline of an essentially cylindrical passage (5.9 cm initial diameter by 

11.1 cm long); thirty of these passages are arranged at equal azimuthal positions 
at the upstream end of a single, annular burner (Fig. 7.1). A 15-deg swirl angle 
is imparted to the air. The fuel flow rate produces an overall equivalence 
ratio of 0.6. 

The Annular Diffuser Deck (ADD) code was used to calculate the air flow 
field within the swirl tube using the assumed inlet air axial and tangential 
(swirl) velocity profiles shown in Fig. 7.2. 

The Particle Tracking (PTRAK) code was used to estimate the spatial dis- 
tribution of fuel vapor sources due to the vaporization of injected fuel. The 
fuel was divided into five classes distinguished by the the angle of injection; 
each class contained twenty percent of the injected fuel. The Sauter mean 
diameter (SMD) for each class was 35 ym. Data for the pressure atomizer used in 
the swirl tube design indicate that the fuel spray is a hollow cone of 80-deg 
included angle with a variation of t 10 deg. For the simulation, the droplets 
were assumed to be injected with angles of 30, 35, 40, 45 and 50— deg relative 
to the centerline of the passage. A fuel injection velocity of 17.6 m/sec 
(corresponding to the specified injector pressure drop of 1.22 atm) was assumed. 
(Because the initial diameter of the droplets is small, the droplets will 
accelerate rapidly to the freestream velocity so that the calculated results are 
insensitive to the choice of initial velocity in these cases.) 

The calculations were performed usinj a 30-deg segment of each swirl tube. 
The fuel droplets acquire a tangential velocity component due to the swirl 
velocity of the air. The PTRAK code accounts for droplets exiting through one 
"boundary" of the segment by allowing identical droplets to enter at the opposite 
boundary of the segment. The radial and azimuthal variations of droplet position 
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for the droplet injected at an angle of 40 deg relative to the axial direction 
are shown in Fig. 7.3 where it can be seen that the fuel droplets are 
to the inner region of the flow. The calculations indicate that over 99 percent 
of the fuel is vaporized at the exit of the swirl tube (Fig. 7.4). 

The Vapor Diffusion (VAPDIF) code was used to calculate the fuel vapor 
concentration profiles throughout the segment of the swirl tube. A turbulent 
Schmidt number <Sc t - u/p 2> t )of unity »as assumed. A contour plot ° f ^uiv.lenc 
ratio at the exit of the swirl tube (and, hence, entrance to the annular burner) 
is shown in Fig. 7.5; it is noted that little diffusion of 

occurred. There is published evidence to suggest that the tubulent Schmidt number 
should be about 1.0 (Ref. 7.1); however, unpublished Corporate-sponsored analytical 
and experimental efforts at UTC suggest that a lower value (0.5) may be more 
appropriate. A lower Schmidt number cotresponds to a higher rate °£ <>iffus . 

In addition, the average value of eddy viscosity calculated by the ADD code is 
lower (bv about a factor of 3.0) than anticipated based upon the use of other 
eddy viscosity models at UTC. The uncertainty in the calculated 
ratio contours due to the uncertainty in the proper choice of turbulent Schmidt 
number and eddy viscosity can be reduced only through the conduct of experiments 
to calibrate properly the model. 

The available evidence indicates that the rate of mixing of the fuel vapor 
and air is probably understated in these calculations. However, if the ealeu 
prof iles'of ^equivalence ratio are at all representative of conditions within a 
swirl tube, then it can be concluded that the design does not Ptoduceaprof lie 
satisfactory for use with the lean, premixed combustion concept. On the other 
hand the design is capable of achieving nearly complete vaporization, a 
result is not expected to be affected by a different choice of mixing parameters 

(u t and Sc t ) . 
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7.2 Scries Staged Premixing Paaaage 


In the series staged configuration, fuel ia injected through sixty fual 
injectors spaced at equal azimuthal locations within a premixing passage that is 
concentric with the axis of the engine (Fig. 7.6). No swirl is imparted to 
either the air or the fuel. The fuel injector is a simple streamlined shape 
(Fig. 7.7) with orifices at four equal azimuthal positions; the fuel is injected 
normal to the local flow direction. The axis of the fuel injector ia parallel 
to the local flow direction. At the cruise condition, the fuel flow rate in the 
premixing passage produces an overall equivalence ratio of 0.55. 


The ADD code was run with the uniform inlet air conditions described above 
and the geometry shown in Fig. 7.6 up to the point where the premixing passage 
just enters the combustor. Four classes of droplets were used in the PTRAK 
analysis. Each class had an initial SMD of 25.14 ym. Each class of droplets 
was injected at an angle nearly normal (80-deg) to the local freestream direction 
but at four equal azimuthal angle locations (corresponding to the four orifice 
positions) around the injector. A 6-deg segment of the annular premixing passage 
was analyzed. The calculations indicate that all of the fuel is vaporized upstream 
of the exit of the passage (Fig. 7.8). 


The fuel concentration profiles were estimated using the VAPDIF code. 

Contour plots of equivalence ratio are shown in Fig. 7.9; as in the case of the 
swirl tube design, it is seen that little mixing of the fuel and air has occurred. 
The turbulent Schmidt number was assumed to be unity and the average eddy 
viscosity calculated by the ADD code was lower than the anticipated level. Thus, 
for reasons similar to those presented in the discussions of the swirl tube design, 
the degree of mixing may be understated in the series staged design. 
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7.4 List of Symbols 
Turbulent mass diffuslvity (cm^/sec) 
Sc t Turbulent Schmidt number u t /pj^ t 

Eddy viscosity (gm/cm-sec) 

P Density (gm/cnr*) 




L TUBE PREMIXING PASSAGE GEOMETRY 




RADIAL POSITION — cm 












PERCENT EVAPORATED 


R82-915362-40 


ORIGINAL PAGE 13 
OF POOR QUALITY 


AXIAL VARIATION IN THE EXTENT OF FUEL EVAPORATED IN THE 
SERIES STAGED PREMIXING PASSAGE 
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3.0 ANALYSIS OF AUTOIGNITION 


8.1 General Approach 

The analysis described in Sections 3.0 through 5.0 provides a method for pre- 
dicting the performance of an LPP design with respect to fuel vaporization and mix- 
ing. This section describes a method for predicting the occurrence of autoignition 
in turbulent two-phase flow within practical premixing section designs. A litera- 
ture review conducted to identify the important characteristics that influence the 
occurrence of autoignition and to aid in the development of autoignition models is 
presented in the Appendix. 

Two autoignition models are described in this section. In either model, it is 
assumed that the progress of the pre-ignition chemical reactions can be related to 
the level or rate of change of concentration of a critical, intermediate chemical 
species. The characteristics of this species and the autoignition criterion are 
specified for each model. As was done in the previous sections of this report, the 
weak interaction assumption is made so that the determination of the occurrence of 
autoignit ion can be carried out in successive steps. First, the viscous flow field 
for the carrier gas (air) is determined by use of the ADD code as described in 
Section 3.0. Second, the liquid fuel droplet heating and evaporation rates are cal- 
culated by the PTRAK code as described in Section 4.0. Third, the fuel vapor con- 
centration distribution is calculated using the fuel vapor source terms and the 
VAPDIF code as presented in Section 5.0. Finally, with the distribution of fuel 
vapor known, the net production and turbulent mixing of critical species in the 
vapor phase can be calculated as described below and an estimate of autoignition 
time in realistic LPP mixing passage can be made. The models described in this 


section permit the calculation of the net production of critical species in a 
two-phase fuel-air mixture in which critical species are produced and consumed both 
throughout the diffusing fuel vapor-air mixture and within the fuel vapor film sur- 
rounding the liquid droplets. This second source of critical species is a function 
of fuel droplet heating and vaporization rates. The basic assumption of weak inter- 
action has beer, applied throughout this formulation of the general problem of pre- 
mixing passage flow analysis. For the autoignition analysis, it is assumed that the 
production of critical species does not affect the fuel vapor concentration. The 
significance of this assumption on the ignition delay times predicted by each of the 
two autoignition criteria is discussed below. An exception to the weak interaction 
has been made in that the tempenture depression due to fuel evaporation is taken 
into account. That is, the enei which is extracted from the carrier gas and which 
is used to heat the fuel droplets id to evaporate the fuel results in a decrease in 
the temperature of the gas. This can significantly affect the pre-ignition chemical 
reaction rates. The amount of heat extracted from the gas is calculated in the PTRAK 
code and must be equal to the change in internal energy of the flow passing through 
a given volume of the flow field. Because the mass flux is also known, the tempera- 
ture change can be calculated. 

Symbols used in this section are defined on pages 114-115. 
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8.2 Basic Equations 

It is sashed that autoi 8 nitio„ is l^tfon!' 

rr.^« is desunatei ° 2) 

is given by the diffusion equation 


-^ 2 C,= S Q +S ( 


( 8 . 2 . 1 ) 


where the convection term is 


DC2-r£^ii$£2. + r^il^2_ + f^_l-^- 

0 1 " L hi •* dy , *- h- ■* dya *■ h 3 ■* dy 3 


( 8 . 2 . 2 ) 


and the diffusion term is 




a 2 C 2 

dy 2 2 


■s c *3 2 


h,h 2 h 3 dy 2 


h,h 2 h 3 


h i^2 

h. 


(8.2.3) 


The terms within the brackets are known minus 

tion 5.0). The source term S g represents the ies in the diffusing fuel 

consumption) per unit time per unit vo ““ the ne t rate of formation per unit 

ir^srs-- «- — - 

.lroplat ■ The source term S g is calculated from 


S9= 


(8.2.4) 


t ion tor critical species are related by 


(8.2.5) 


e- 

. • j Y /j*- nvpr the t irne At in which the 
S ovlr the volume of vapor film surround- 

Utho droplet, and summins over all droplet classes. 
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v*. 


,r«+b 
d 4nr* 


,t+At 


dt dr} 


( 8 . 2 . 6 ) 


Eq. (8.2.6) is evaluated at the local temperature and fuel concentration (partial 
pressure) of the vapor surrounding each droplet. The radial variations in tempera- 
ture and fuel vapor partial pressure are determined in accordance with the model 
presented in Section 4.0 which is based upon the model described in Ref. 8.1. The 
radial variation in film temperature is given by 




(8.2.7) 


It can be shown that the Lewis number within the fuel vapo:: film is generally not 
equal to unity so that the thicknesses of the heat and mass transfer boundary 
layers surrounding the droplet differ. However, the heat transfer and mass trans- 
fer coefficients used in the model developed in Ref. 8.1 are determined from corre- 
lations that use droplet diameter as the length scale; this convention has the 
practical effect of defining the heat and mass transfer boundary layer thicknesses 
as equal. Essentially, the model presented in Ref. 8.1 is used to derive a func- 
tional form for the vaporization rate and then a mass transfer coefficient corre- 
lation is used to obtain the constant of proportionality. Consistent with this 
approach, the functional form for the radial variation of fuel vapor partial pres- 
sure is derived assuming that the heat and mass transfer boundary layer thicknesses 
are equal. 


P f= p o-t p o- p f,s)exp[A(r-r # )] 


( 8 . 2 . 8 ) 


where 


A * — ■ — in( — ~ — ) 

s V p o-Pf,s 


(8.2.9) 


Then the molar concentration of fuel is: 


*i sp f/(RoV 


( 8 . 2 . 10 ) 


The decrease in temperature of the carrier gas (air) due to the heat extracted 
by droplet heating and vaporization is calculated from the energy equation. If it 
is assumed that diffusion and crosswise convection can be neglected, then 


.l-SL.- 

J <iy, 




( 8 . 2 . 11 ) 


where the sink term on the right-hand side is obtained by integrating the droplet 
heat transfer rate ever the time the droplet resides in the volume dV and then summing 
these integrals for all droplet classes. 
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8.3 Autoignition Models 


The first autoignition model described in this section (designated Model I) is 
based on a kinetics scheme proposed by Hautman, et al . (Ref. 8.2) who examined data 
for the autoignition of several paraffin fuels and determined that during the pre- 
ignition period (1) the fuel molecule decomposes into ethene and other chemical 
species and (2) the maximum ethene concentration corresponds in time to the time 
at which a sudden rise in gas temperature occurs. Thus, the time of maximum ethene 
concentration is defined by the authors as the autoignition time. It is assumed in 
the model described in Ref. 8.2 that all of the fuel is in the vapor phase, that 
the amount of fuel is specified as an initial condition, and that all chemical reac- 
tions occur in the vapor phase at the molecular level. It is assumed herein that 
this model can be generalized to two-phase flows consisting of mixtures of liquid 
fuel droolets and fuel vapor convected by a carrier gas such as air. The rate of 
change of the molar concentration of ethene is given by 


dt 2 R F X 0 2 ''I "'2 -'B-Og -I -2 (8.3.1) 

(where X-^ and X 2 represent the concentrations of fuel and critical species, respectively) 
and the reaction rate constants are of the form 


ap 


- R„X 


qb /9b ye 


R f = a f exp(-E F /R*T) 


(8.3.2) 


R b = Ag exp(-E a /R*T) 


(8.3.3) 


Although Eq. (8.3.1) has the formal structure consistent with forward and reverse 
reactions, it could represent the net production and depletion of ethene by 
distinct processes. In terms of mass fraction, Eq . (8,3.1) becomes 


S - — p X^C^C^-R y aB c 


(8.3.4) 


where 


_ r f 

Rf '' 


(8.3.5) 


/f/ 2 /0^ B+yB R B 

R ®' 


(8.3.6) 


In the model described in Ref. 8.2, the rate of depletion of fuel is calculated 
according to: 
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* -R. 




(8.3.7) 


Eqs . (8.3.1) and (8.3.7) are solved simultaneously. Since 8g is negative, then the 
ethene concentration will (eventually) reach a maximum and decrease thereafter. The 
time of maximum ethene concentration (as indicated by dX 2 /dt equal to zero) is 
designated as the autoignition time. 


Due to the weak interaction assumption, no fuel depletion chemical reaction is 
included in the present analysis. Thus, the concentration of ethene will approach 
asymptotically a maximum value and an autoignition criterion different from that 
used in Ref. 8.2 must be developed. In the present analysis, autoignition is assumed 
to occur when the time rate of change of ethene due to chemical reaction becomes 
less than some specified fraction of the local ethene concentration; i.e.. 


dx 2 /dt 

< ui (8.3.8) 

x 2 

It has been estimated by the present authors using the data presented in Ref. 8.2 
that a) is of the order of 5000 sec - ^. A refined estimate of u> would require calibra- 
tion of the criterion used herein with the numerical results reported in Ref. 8.2 

It should be noted that the model described in Ref. 8.2 is based upon data ob- 
tained using a shock tube operated at stagnation temperatures typically in excess 
of 1200K whereas the present analysis is intended for use in the analysis of pre- 
mixing passages (stagnation temperatures between 400 and 900K) . The constants used 
in the model described in Ref. 8.2 were obtained by calibrating the model with these 
shock tube data. It is well-known that different rate-controlling processes may 
be important at high and low temperatures. The present authors estimate that the 
use of the model constants presented in Ref. 8.2, when applying the model to the 
analysis of flows within premixing passages, will produce unrealistically large 
values of ignition delay time. Clearly, calibration of this model for the lower 
temperature (and higher pressure) levels characteristic of flows within premixing 
passages is required and therefore, there is no basis for selecting a value of the 
autoignition criterion (w) for these conditions. 

In the second model for autoignition (Model II), it is assumed that autoigni- 
tion occurs when the concentration of an unknown, critical species reaches a criti- 
cal value. This approach is essentially the same as that described in Ref. 8.3. 

Model constants are obtained by comparing the calculated autoignition time with ex- 
perimental data. Thus, if X 2 represents molar concentration of critical species 
and if X 2 C is its critical value, then ignition occurs when 


x 2 sx 2c 


(8.3.9) 
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or equivalently. 

X 2 

X 2C *' (8.3.10) 

It is assumed th.it the normalized rate ot product ion of critical species is given 
hv an expression of the form: 


d(x 2 /x 2c ) 

dt 


- A F e 


- Ef/R * T * ( x 02 *V 


(8. J. 11) 


and tiuit tin* production of X> has nop l i y % i b 1 e effect on t ho concent rat ion of fuel 
.uni ox i d i /.o r . 


In ordet t o nst' 1 q . ( S . \ . I l) in tin' diffusion equation, Kq. (S.d.l), it is 

mvi'ss.irv either ( I ) to know tin' value ot X> c or (12) to express both the source 
terms and tin 4 ditiusiou tMpwit ion in tonus of normalised conceit t r a t i oils . because 
in this model t ho charac t or i s t i os o t tho critical species arc unknown, tho socotul 
approach is takon. Tho source terms >\iven hv Kqs . (8.2.4) and (8.2 6) can bo ox- 
p i e s se d a s : 


S 


q 


= X 2C 



18 . 1 . 12 ) 


and 


s cO 


x 2c f< 


i[d(x 2 /x 2c )/dt 


(s.i.io 


and tht' d illusion equation can ho wr i t t on as: 


DIC 2 /C£c> 

Dt 


V (C^/C^) 


: 2C 




v 8 . 1.14) 


\ 




ilii'li , us l Ur. Ii| . is . > I . 


0(Cp'C ?c ) 

Dt 


■^<C 2 /C 2c ) 


AT'VV 


(8 . I . IS) 


ft ) •; * “ \ i dt ut t i on examination ot ftp; ( S . 2 . 4 ) , ( 8 . 2 . h) . and (8.1. IS) that tht' 
i i i'i 1 1 i\l irlntive ria-o; ! i m t i on distribution is i tulcpetulen t ot tin* molecular 

u : ■ i ’ll! , * ! thfiiitii.il * pOi it**; . 

!h. i I'lisi.m! •; ap t .ai tuc. in lq. (8. I It) are obtained hv lirst deriving an 
rxj'irea.'ii tot tin' ten it ion delav tine t cr a uniform fuel air mixture in a constant 
t it - ne i a t ut e , constant ptessute t * ow and then comparing the tesult t o the appropriate 
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ignition delay time correlation for this aimple flow situation. The autoignition 
time t is defined such that 


or 


and therefore 



r T d(X 2 /X 2c ) 

l di — 




(8.3.17) 


T s 


A -! e Ef/R* T 


(8.3.18) 




t 

i 


! 


As an example of the method used to determine the constants in Eq. (8.3.11), 
consider the ignition delay time data for various fuels obtained by Spadaccini and 
Te Velde (Ref. 8.4) and recently correlated in terms of mixture temperature and 
equivalence ration (Ref. 8.5). A suggested correlation assuming second-order pres- 
sure dependence is, for example only: 

xe E/RT 

T= p2<£ (8.3.19) 

The temperature appearing in Eq. (8.3.19) is the mixture temperature corresponding 
to the equivalence ratio, 4> . The equivalence ratio is related to the molar concen- 
trations of fuel and oxygen by 


so that 


and 


<*> = 


Xf/X 02 


(8.3.20) 


*f = X 0 2 ^ n ST 


(8.3.21) 


*o 2 = yo 2 


P 

Rq T 


(8.3.22) 


Thus, Eq. (8.3.19) becomes 


*o 2 "st 

2 

Ap 0 2 


£./*• T 

e 

P 2 <f>' + * T 2 ** 


(8.3.23) 
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Then, comparing Kqs . (8.3.19) and (8.3.23), the constant used in the model are 


Yt K 


Er = F, 


b = 0 


A 
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8.4 Muter ical Methods 


The source term S d given by Eq. (8.2.6) and the energy sink tern in Eq. (8.2.11) 
ere calculated in the PTRAK code using a predictor-corrector method described in 
Section 4.11. This method ie second-order accurate over the time step dt. The 
volume integration in Eq. (8.2.6) is obtained using the trapesoidal rule which is 
also second-order accurate. 


The solution of the non-linear diffusion equation, Eq. (8.2.1), is obtained 
using the numerical methods described in Section 5.3 and is included as part of the 
VAPDIF code. The source term S g requires special treatment when Model 1 is used 
since Eq. (8.3.4) is nonlinear in C 2 . The authors of the model described in Ref. 8.2 
recommend an integration step-sise on the order of one microsecond whereas the typi- 
cal integration step-sise (in terms of flow residence time) used by the VAPDIF code 
is at least one order of magnitude larger. Presently, the source term S g is deter- 
mined from 


s 

• 2 At 


(8.4.1) 


where 


AX2_ = J_ r dXe_ 

At At J m dt 


(8.4.2) 


such that the integration step-sise is approximately 1 nsec in the volume whose 
residence time is At. 

Since the source term is not a function of C 2 (or equivalently, X 2 ) for Model 
II, it is not necessary to calculate S g using Eqs . (8.4.1) and (8.4.2). Instead, 
a single calculation is performed using Eqs. (8.3.11) and the source term is then 
determined using Eq. (8.2.4). For moderate to large numbers of computational grid 
points, the use of Model II can result in a significant savings in computational 

tisie . 
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8.5 Results and Discussion 


Autoignition Model 1 was used to analyze the flow within the Series Staged 
Premixing Passage described in Section 7.2. The constants used in this model were 
obtained from Ref. 8.2 and are presented in Table 8.1. Figure 8.1 shows the cal- 
culated axial variations of ethene mass fraction along several coordinate lines, 
which in the present approach approximate the actual streamlines. These stream- 
lines are located at approximately 20, 40, 60 and 80 percent of the distance from 
the inner to the outer wal’ of the passage in the plane of symmetry. The results 
presented in Fig. 8.1 show a rapid increase in ethene mass fraction in the region 
where the fuel droplets are evaporating. The ethene mass fraction increases at a 
less rapid rate as droplet vaporization is completed and then appears to increase 
more rapidly near the duct exit where only the gas-phase reaction occurs. These 
calculated results were obtained for two-phase flow and include the effects of 
droplet evaporation, turbulent diffusion, and chemical reaction both in the gas 
phase and in the film surrounding the liquid fuel droplets. 

The autoignition criterion proposed for Model I, Eq. (8.3.8), states that 
autoignition occurs when the rate of change of ethene concentration due solely to 
chemical reaction becomes less than some fraction of the ethene concentration. 

The axial variation of the minimum rate of change in ethene concentration due to 
chemical reaction relative to the local ethene concentration is shown in Fig. 8.2. 

For a completely premixed, preveporized fuel air mixture at the point of injection, 
as assumed in Ref. 8.2 and as assumed in deriving the autoignition criterion, the 
autoignition criterion should decrease monotonically . However, the mixture within 
the Series Staged duct is neither completely prevaporized nor premixed at the point 
of injection so that the behavior of the autoignition parameter shown in Fig. 8.2 
is quite different. Furthermore the absolute level of this parameter is much less 
than that suggested by examining its behavior for the relatively high temperature 
flow situations described in Kef. 8.2 Therefore, without further calibration it 
is not possible at the present time to use Model 1 and autoignition criterion, 

Eq. (8.3.8), to determine the likelihood of autoignition in the Series Staged passage. 


Autoignition Mod-1 11 was then used to analyze the flow within the Series Staged 
Promixing Passage (described in Section 7.2). The constants used in this model 
were obtained by applying the procedure described in Section 8.3 to an unpublished 
c- relation of the ignition delay data for Jot— A fuel presented in Ref. 8.4. 


r = 


8 59xiO- ,2 e 366,4/R#r 

p 2 4 > 


(8.5.1) 


The constants for Model II are presented in Table 8.2. Calculations made using 
this correlation indicate that autoignition within the Series Staged Premixing 
Passage does not occur for the assumed flow conditions if the flow therein is uni- 
form. The contours of relative mass fraction, calculated using Model II and shown 
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in Fig. 8.3, indicate that autoignition may occur near the exit of the duct in the 
relatively fuel-rich interior of the flow (see Fig. 7.9). Since the calculated 
profiles of equivalence ratio and mixture temperature are non-uniform, a direct 
comparison of calculated autoignition times made using Model II and the correlation 
based upon uniform flow conditions (Eq. 8.5.1) is not possible at present without 
additional calibration of the model. 



112 


R82-915362-40 


8.6 References 

8.1 El Wakil , M. M. , 0. A. Uyehara and ?. S, Myers: A Theoretical Investigation 

of the Heating Up Period of Injected Fuel Droplets Vaporizing in Air, NACA 
TN 3179, May 1954. 

8.2 Hautman, D. J., F. L. Dryer, K. P, Shug and I. Glassman: A Multiple Step 

Overall Kinetic Mechanism for the Oxidation of Hydrocarbons, Combustion Science 
and Technology, Vol. 25, 1981, pp. 219-235. 

8.3 Faeth, G. M. and D. R. Olson: The Ignition of Hydrocarbon Fuel Droplets in 

Air, SAE Paper No. 68045, May 1968. 

8.4 Spadaccini, L. J. and J. A. TeVelde: Autoignition Characteristics of Aircraft- 

Type Fuels, NASA Contractor Report CR-159886, June 1980. 

8.5 Spadaccini, L. J. and J. A. TeVelde: Autoignition Characteristics of Aircraft- 

Type Fuels, Combustion and Flame, 1982 (to be published). 



113 


R82-915362-40 






A 


A f , a b 


b 

C P 

C 

Ep i Eg « E 
h 


h\ ,h 2 ,h 3 

K 


M 


n 

n st 

P 


Rp» R B* R F» R B 


R o 

R* 

S C 

s H ,s 

d g 

t 

T 

Ul,U 2 ,U 3 

V 

X 


8.7 List of Symbols 

Constant (1/cm) 

Reaction rate constants 

Film thickness = droplet radius (cm) 

Specific heat at constant pressure (cal/gm/K) 

Mass fraction 

Activation energy (cal/mole) 

Heat transfer coefficient (cal/cm^/K/sec) 

Metric scale coefficients ( 1/V, 1/V, r, see Section 3.7) 
Constant 

Molecular weight 

Droplet number density (1/sec) 

Stoichiometric oxidizer to fuel mole ratio (dimensionless) 
Pressure (atm) 

Reaction rate constants 

Universal gas constant (82.0575 cm^-atm/mole/K) 

Universal gas constant (1.98717 cal/mole/K) 

Schmidt number (dimensionless) 

Source terms for critical species (gm/cm^/sec) 

Time (sec) 

Temperature (K) 

Velocity components (m/sec) 

Metric scale coefficient (inverse of potential flow velocity) 

<■* 

Molar concentration (moles/cnr ) 
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Coordinates (dimensionless) 

Mole fraction of oxygen (dimensionless) 

Blowing parameter 

Reaction rate exponents (dimensionless) 

Residence time in volume element (sec) 

Effective turbulent viscosity (gm/cm/sec) 

Carbon number of fuel (dimensionless) 

Equivalence ratio (dimensionless) 

Density (gm/cm^) 

Residence time (sec ) 

Autoignition criterion (sec~*) 

Air 

Critical value 
Fuel 

"Reverse" reaction - depletion of critical species 
"Forward" Reaction - production of critical species 
Oxygen 

Droplet surface 
Fuel 

Critical spec ies 
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TABLE 8.1 

Autoignition Model I Constants Used in the 
Analysis of the Series Staged Premixing Passage 


n 

12 



a f 

1.0 x 10 17 * 32 

b B 

1.0 X 10 

e f 

49600 


50000 

a 

1.07 

a B 

1.18 

6 

.5 

6 B 

-.37 

Y 

.4 

7 b 

.90 


TABLE 8.2 

Autoignition Model II Constants Used in the 
Analysis of the Series Staged Premixing Passage 


n 

A 

E 

a 

3 

<S 


t 

l 


12 

3.05 x 10 17 
36614 
1.0 
1.0 
0.0 
2.0 
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9.0 CONCLUDING REMARKS 


1. The analytical approach for modeling the operational characteristics of 
premixing-prevaporizing fuel-air mixing passages developed herein (i.e., 
the successive application of the ADD, PTRAK, and VAPDIF computer codes) 
produces results which are physically realistic. Insufficient data exists 
for calibrating the fuel droplet evaporation model (PTRAK) and the vapor 
diffusion model (VAPDIF). Predictions of the evaporation model, however, 
are consistent with the available data. Results of calculations performed 
by the diffusion model indicate that the turbulence levels and diffusion 
rates are lower than levels anticipated on the basis of general experience. 

Data on turbulence levels and turbulent dif fusivities are required to cali- 
brate the turbulence models. 

2. The model sensitivity study indicates that the more important parameters for 
determining the fuel droplet evaporation rate are initial droplet size, fuel 
volatility, and ambient temperature, while the less important parameters are 
initial fuel temperature and initial droplet velocity. 

3. The model was used to analyze two designs of premixing-prevaporizing fuel-air 
passages. The model predicted that both designs produce essentially complete 
vaporization. For the levels of turbulence calculated by the model, it was 
determined that neither design produces a uniform vapor ft el-air ratio profile 
at the exit of the passage. 

4. Two autoignition models have been described and these models have been incor- 
porated into the analysis. One model embodies the framework required for the 
incorporation of multi-step chemical kinetic reaction mechanisms. The other 
is based on application of a simple global expression for describing the pro- 
gress of the autoignition process. Both models were applied to the analysis 
of a premixed-prevaporized fuel-air mixing rassage; however, only a qualita- 
tive assessment of the computed results was possible since further calibration 
of these models is required before a quantitative assessment can be made. Pre- 
sently, a single rate expression is employed in the first model; with further 
development of the VAPDIF code, simplified multi-step mechanisms for hydrocarbon 
fuel oxidation appearing in the literature could be incorporated. In addition, 
further calibration and refinement of the autoignition models used in the 
present calculations are needed which account for the depletion of fuel and for 
the influence of turbulence on pre-ignition reaction rates. 
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APPENLIX A - LITERATURE SURVEY 


LITERATURE SURVEY CONDUCTED IN SUPPORT OF 
THE DEVELOPMENT OF A MODEL FOR AUTOIGNITION 


INTRODUCTION 

A literature survey was conducted to provide information in support of the formula- 
tion and calibration of an analytic model capable of predicting the onset of autoxgnr- 
tion within the premixing passage of lean, premixing, prevapor is ing (LPP) combustors. 
These models have been incorporated into the computer programs used to analyte the LPP 
combustor preroixing passage flows. These computer codes are described in the previous 
sections. The results of the literature survey are presented herein. 

Symbols used In this section are defined on page 136. 

Background 

The application of the LPP concept presents several design problems a major 
problem being prevention of autoignition. It ia well known that autoignition time 
(ignition delay time) decreases with increasing pressure and gas temperature. There- 
fore. the severity of the autoignition problem increases as the design pressure ratio 
of modern gas turbine engines increase. For most proposed designs, autoignition in 
the premixing passage must be prevented because the uncooled engine hardware in this 
duct can undergo catastrophic failure aa a result of the sudden large heat release 
following autoignition. Clearly, the exploitation of the LPP concept requires that 
the residence time within the fuel preparation passage be long enough to achieve 
essentially complete vaporisation and mixing of the fuel and yet short enough to pre- 
clude the occurrence of autoignition. The analysis described in Sections 3 thorugh 5 
( l *>80) provides a method for predicting the performance of an LPP design with respect 
to fuel vaporisation and mixing. It is desired to incorporate within the analysis a 
model of the autoignition process such that it will be possible to determine whether 
autoignition occurs in a given fuel preparation duct. This model must be sufficiently 
comprehensive to address all of the processes which could affect the progress of the 
pre ignition reaction in the turbulent, two phase flow field characteristic of prac- 
tical premixing section designs. 

The first objective of this effort was to review the literature from the point 
of view of identifying the characteristics of the flow which can influence the 
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occurrence of autoignition. The following factors were projected to be of signifi- 
cance: 

mixture pressure; 
mixture temperature; 

degree of fuel vaporization and mixture ratio; 
droplet initial temperature; 
droplet size; 

mixture turbulence characteristics including the intensity, scale and 
distribution of turbulence; 
fuel chemical composition. 

Attention was given to examining the literature from the standpoint of defining the 
appropriate functional relationships between autoignition time and these factors. 
Attention was also given to literature which might point to factors of significance 
other than those listed above. 

When conducting the survey, primary consideration was given to papers dealing 
specifically with autoignition or ignition delay phenomena. However, related infor- 
mation may prove useful in developing this model; therefore, information on the 
autoignition temperature, the mixture below which autoignition is not possible, and 
on the minimum energy which an external source (e.g., a spark) must apply in order 
to ignite a fuel-air mixture is cited where such information is believed to be of 
significance . 

The experimental procedures employed in various efforts to obtain autoignition 
data were not directly of interest in this survey; the reader interested in these 
techniques is referred to Spadaccini and TeVelde (1980). However, as emphasized by 
Spadaccini and TeVelde, the interpretation of reported data does, in fact, require 
an appreciation of the influence of the experimental technique used upon the data 
obtained. 

The review of the experimental literature on autoignition is given in the fol- 
lowing section of this report. A second objective of this task was to review previ- 
ously published models of the sutoignition process. The results of that phase of 
the effort are given in the final section of this appendix. 

Procedures 

This survey was conducted by reviewing the contents of the United Technologies 
Corporation Library. In addition, computer-assisted searches were conducted of the 
National Technical Information Service, Science Citation Index, and Current Contents. 
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REVIEW OF EXPERIMENTAL DATA 


1. Effect of Mixture Pressure 


It is reasonable to expect that the rates of gas phase chemical reactions vary 
directly with the concentration of the reactants. For pre-ignition reactions of 
air and fuel in the range of mixture ratios of interest, it follows that reaction 
rates should vary directly with pressure. Because ignition delay time is related 
inversely to the pre-ignition reaction rates, ignition delay time will vary inversely 
with mixture pressure. Likewise, autoignition temperatures and minimum ignition 
energy should vary inversely with mixture pressure. 

Many ignition delay time experiments have involved injection of liquid fuels 
into hot gas streams; therefore, to interpret the results in terms of the pressure 
effect on gas phase reactions, consideration must be given to how pressure variations 
could influence the results other than through the effect on gas— phase reaction rates. 
In particular, the extent of fuel vaporization of the liquid spray is a function of 
mixture pressure. As pressure increases, droplet vaporization rates are influenced 
by two factors: 1) enhanced heat transfer from the gas phase to the liquid phase, 

and 2) decreased heat transfer due to a reduced temperature differential when the 
droplet reaches its higher equilibrium temperature. At the same time, the drag on 
the droplet due to the relative motion between the droplets and the gas increases 
because of the increased dynamic pressure. Assuming that the droplet is injected 
with a velocity less than that of the air stream, then the residence time of the 
droplet within a fuel preparation passage of fixed length decreases with increasing 
pressure. The amount of fuel vaporized is the integrated product of vaporization 
rate and residence time. Calculations performed at UTRC for representative condi- 
tions indicate that the amount of fuel vaporized within a fixed duct length changes 
only slightly with changing pressure [see, for example, Anderson, et al. (1980)]. 

It is therefore reasonable to assume that the pressure effects reported in the lit- 
erature for experiments conducted with both liquid and gaseous fuels can be attributed 
primarily to the effects of pressure on gas-phase kinetics. 

Correlations of the effects of pressure on ignition delay are usually expressed 
in the form: 

T = 1/P n (1) 

A summary of the correlations found in the literature for fuel injected into flowing 
air streams is presented in Fig. A.l. The number in parentheses next to each data 
correlation is the value of the exponent, n. Spadaccinl and Tevelde (1980) provided 
correlations of their data using values of the exponent n equal to 1 and 2; both 
sets of correlations are shown in Fig. A.l. As can be seen, the body of data available 
in the literature indicates a suitable value for n lies in the range from 1 to 2; a 
more precise generalization cannot be made. 


R82-915362-40 


Numerous experiments on related phenomena supported the findings reported for 
the flowing air stream experiments. For example, Halstead, et al. (1977) report ig- 
nition delay data for various liquid fuels in reciprocating engines. These data 
indicate that ignition delay time varies inversely with mixture density to the first 
or second power. Since both pressure and pre-ignition delay temperature vary during 
the compression stroke, the effect of pressure on the ignition time data reported 
from such experiments is not easily determined. However, the variation of ignition 
delay time with density implies a variation not inconsistent with the pressure be- 
havior presented in Fig. A.l. Data for the variation of autoignition temperature re- 
ported by Ingebo (1974) for fuel-air mixtures in a combustor environment and for 
minimum ignition energy reported by Satcunanathan (1971) for individual droplets in 
an oxidizing environment also show an inverse correlation with mixture pressure. 


2. Effect of Mixture Temperature 

It is generally accepted that gas phase chemical reaction rate behavior can be 
represented using reaction rate constants having the form: 

i * -E/RT 

k * Ae (2) 

This form of expression reflects the fact that as mixture temperature increases, the 
kinetic energy of the reactant molecules increases, intermolecular collisions become 
more energetic, and the probability of reaction increases. Indeed, it has been the 
practice by most experimenters to use expressions of the form: 


to correlate experimental results. It is important to recognize that the value of the 
activation energy, E, reported by the experimenters is in reality a correlation co- 
efficient and not the activation energy for a specific reaction. Thus, if the ex- 
perimental data is to be used to specify the constants appearing in a global or 
single-step rate expression, the form of the rate expression must reflect the form of 
the expression used to correlate the experimental data. Specifically, if an activa- 
tion energy is obtained from a 9et of experiments in which the data is correlated 
using the form 

t * ap e (4) 

then the reported activation energy should only be used in an expression which shows 
a second order pressure dependence. For example, Spadaccini and TeVelde (1980) report 
values of activation energy obtained from regression analysis employing either first 
or second order pressure dependence and show that the activation energy values ob- 
tained assuming first order dependence are approximately 15 percent lower than those 
obtained assuming second order dependence. 


It is also important to recognize that the experimentally-determined values of 
activation energy in general reflect phenomena other than gas-phase chemical reactions 
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v* *° * h V mperf *? t n * tur * of th# « x P«riment* . Obviously, for those experiments in 
which liquid fuel is injected into a flowing air straaa, increased mixture teapera- 
tura affect, droplet y.poriaation ratea and therefore the concentration profilea be- 
tween the point of injection and the point of ignition will differ for different 
air temperatures. 


Data from a representative group of experiments is presented in Fig. A.2. The 
slope of each correlation curve is a measure of the activation energy; the calculated 
value of activation energy is given by the value in parentheses adjacent to each 
curve. Note that the data group corresponds to those experiments which were corre- 
lated assuming a pressure exponent, n, of 1.0. A wide variation in magnitude of 
activation energy is observed and is probably wider than the variation due to fuel 
chemistry alone. Differences in experimental technique which result in different 
rates of fuel-air mixing and different rates of spray variation certainly influence 
the experimental results. An example is the work of Myasaka and Mizutani (1975) 
where the apparent activation energy is reported as a function of fuel injection 
rate. Because experimenters are unable to separate these effects from the gas-phase 
kinetic effects, the amount of experimental bias is unknown in all cases. The be- 
havior of the reported data indicate that a dependence of ignition delay time on 
temperature of the form given by Eq. (3) is suitable for modeling, but the appro- 
priate activation energy can only be determined by conducting calibration experiments 

using an apparatus similar in configuration to the fuel preparation system of inte- 
rest. 



An important finding reported by numerous researchers is that although expres- I 

sions of the form of Eq. (3) are useful for correlating data, the range of tempera- ■ I 

tures over which a given correlation applies is limited due to changes in the rate- 
controlling reaction steps as temperature changes. For example, Yashizawa, et al. 

(1978) present data and a reaction rate expression that shows that different values 
of activation energy apply above and below approximately 1200K for butane, hexane 
and octane. Myers and Bart le (1969) show similar results for propane. Brokaw (1965) 
has examined the reaction system responsible for the exponential growth of OH radi- 
cals during the ignition delay period and provides illustrations of how the changing i 

dominant reactions control delay in both the long (millisecond to second) and short 
(less than millisecond) ignition delay regions. The important point is that corre- 
lations or reaction rate mechanisms developed from high temperature experiments can- 
not be used at low temperature conditions without the risk of serious errors. 

3. Effects of Degree of Vaporization and Mixture Ratio ’ 

It is clear that the increasing extent of vaporization of a fuel spray as time j 

increases will affect the rate of progress of the ignition delay reactions through I 

the effect of increasing reactant concentrations on reaction rates. These concen- ! 1 

t rat ion effects are treated as pressure effects where the pressure represents the ! \ 

partial pressure of the reactants. Additionally, as the spray Vaporizes, the fuel . ; 

to oxidizer ratio increases; the mixture ratio may or may not have an effect on re- 
action rates. Also, within a vaporizing spray, two phases are present. No evidence ! 

exists in the literature, however, which would indicate that heterogeneous reaction i 

i : 

i 
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phenomena are of significance. As a result, the only new effect which the degree of 
vaporization may introduce is the effect of mixture ratio. 

With respect to mixture ratio effects, Yashizawa, et al. (1978) obtained igni- 
tion delay time data in gaseous mixtures of n-butane, n-hexane, and n-octane for 
equivalence ratios between 0.05 and 1.0 and determined that ignition delay time did 
not vary with equivalence ratio. Correlations of ignition delay time with equiva- 
lence ratio have also been reported by Marek, et al. (1977) and Halstead, et al. 

(1977) and Spadaccini and TeVelde (1980) for liquid fuel injection. However, as 
noted in these works, the temperature of the fuel vapor-air mixture decreases within 
the vicinity of the fuel droplets due to fuel vaporization; and it is apparently the 
effect of this mixture temperature decrease on ignition delay time which is expressed 
by correlations of ignition delay time with equivalence ratio. Spadaccini and TeVelde 
(1980) presented ignition delay time data both as a function of inlet air temperature 
and as a function of the mixture temperature assuming complete mixing and vaporization; the 
data for each equivalence ratio correlate with the mixture temperature corresponding to 
that equivalence ratio. On the basis of the limited data available, it is concluded 
that, for the range of mixture ratios of interest in LPP systems, mixture ratio does 
not have a large effect on the autoignition process. 

4. Effect of Fuel " itial Temperature 

If the liquid fuel is heated prior to injection, then the fuel will vaporize 
more rapidly than would unheated fuel. Also, the mixture temperature decrease 
associated with a given degree of fuel vaporization will be less for the heated 
fuel than for the cold fuel. Based on the remarks in the preceding sections, igni- 
tion delay time would be expected to decrease with increasing fuel injection tem- 
perature. The data of Spadaccini and TeVelde (1980) for Jet— A fuel preheated to 
400K prior to injection show little effect of initial fuel temperature on ignition 
delay time. The lean, premixing prevaporizing concept requires that essentially all 
of the fuel vaporize within the premixing passage. Achievement of a high degree of 
fuel vaporization within is dependent on the atomization of the fuel into small 
droplets (d < 100 urn). Small droplets will be heated rapidly by the air stream so 
that the effect of injection temperature should be minimal. The estimated Sauter 
nean diameter of the Jet-A fuel injected in the tests described by Spadaccini and 
TeVelde (1980) was approximately 30 pm; these droplets are small enough to be heated 
rapidly by the air at conditions representative of those that exist within the LPP 
passage and therefore the reported results are reasonable. 

The data reported by El Wakil and Abdou (1966) shows about a two-fold decrease 
in ignition delay times as the fuel injection temperature i 9 increased from 330K to 
380K. However, these data were obtained using droplets with initial diameters of 
1650 i.m; droplets of this size are not characteristic of LPP system sprays. Pre- 
heating of the fuel in the case of large droplets helps to overcome the large ther- 
mal inertia of the droplets. 

5. Effect of Droplet Size 

It is reported in the literature that the presence of droplets in a fuel-oxi- 
dizer tream can have a significant effect on the onset of ignition. That is, 
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ignition may occur more readily in a stream containing droplets than a stream of the 
same overall fuel air ratio in which the fuel is completely vaporised. Wood and 
Rosser (1969) have shown that ignition within a system produced by a droplet falling 
through a hot gaseous oxidizer can occur either in the wake of the droplet or in the 
fuel-air mixture envelope surrounding the droplets depending on conditions. With all 
except droplet size remaining the same, ignition delay times remained the same as 
droplet size was decreased until a critical size was reached, then ignition delay 
times increased. It was observed that this increase was associated with a shift in 
the point of ignition from the droplet fuel— air mixture envelope to the droplet wake. 
These experiments indicate that the presence of droplets can have a controlling 
effect and that effect is dependent on the droplet size. The droplet diameters used 
in the experiments ranged in size from 100 to 300 microns. It was concluded that the 
smaller droplets evaporate before ignition occurred in the vapor-fuel/oxidizer mix- 
ture in the droplet wakes. For the larger droplets, pre-ignition reactions progressed 
to the point of ignition in the envelope surrounding liquid fuel droplets; both the 
rate of progress of the pre-ignition reactions and the rate of diffusion of inter- 
mediate species into the environment played a role in this process. 

Other experiments in which the effect of droplet size on ignition delay were 
investigated include the work of El Wakil and Abdou (1966) who determined experimen- 
tally that ignition delay time can be correlated as 


x ~ D. 


for single droplets in the size range from 1200 to 1800 microns for several paraffins 
(n-octane, n-decane, n-dodecane, n-tetradecane , and n-hexadecane) suspended in a hot 
air stream; in these experiments the temperature of the air was varied between 1000 
and 1250K. The authors concluded that the observed dependence on droplet diameter 
could be explained on the basis of consideration of the physical (droplet heating, 
diffusion) and chemical processes occurring in the envelope surrounding the droplet. 
An analysis based on a spherically symmetric field is presented. These authors 
also observed that ignition could occur either in the droplet envelope or in the 
droplet wake; at higher droplet-air velocities, burning occurred in the droplet wake. 
Kadota, et al. (1976) reported no correlation of ignition delay time with droplet 
diameter for single droplets of n-heptane which are plunged into a quiescent furnace. 
The contrasting result obtained by these two experiments with respect to the role 
of the size of the droplet on ignition delay is probably not of significance because 
of the imprecision associated with the small droplet size range studied and the 
difference in experimental conditions. These experiments were carried out with 
isolated droplets in an infinite oxidizer environment. A definitive experiment in 
which autoignition times have been measured for systems with and without presence of 
Jroplots at overall fuel-oxidizer ratios in the range of interest for gas turbine 
combustion systems has not been conducted. 


Effect of Turbulence 


No data were located which show the effect of turbulence on ignition delay time. 
There have been studies which show that increasing turbulence level decreases the 


128 



R82-915362-40 

ease with which ignition is achieved [e.g., Ballal and Lefebvre (1979)] or increases 
the tendency for a flame to blowout [Radhadkrishnan (1979)]. In the case of igni- 
tion, the increased turbulence causes the local rate of dissipation of thermal energy 
to increase relative to the local rate of energy release due to reaction. When dis- 
sipation rate exceeds the heat release rate, propagation of the reaction front, 
which is required for ignition, cannot occur. In the case of blowout of a premixed 
flame, the role of turbulence is to exchange cold reactants with hot combustion pro- 
ducts. Because reaction rates are exponentially dependent on temperature and de- 
pendent on reactant concentration to a lesser degree, the increased rate of exchange 
of cold reactants with hot products always causes a decrease in heat release in the 
flame stabilization wake or pilot region; hence increased turbulence is not favorable 
to blowout limits. In both the ignition and blowout cases, the essential feature 
of the flow is the non-uniformity of the temperature field which the turbulence tends 
to dissipate. In the case of the LPP system, the temperature non-uniformities will 
be a result of non-uniform work performed on the gas during the compression process 
and the temperature depression associated with the injection of cold fuel into the 
hot air. Gradients in the temperature field due to these effects should be small 
and hence the influence of turbulence on redistributing the thermal energy will be 
small. On the other hand, the concentration gradients resulting from imperfect pre- 
mixing can be large both on the scale of the duct and on the scale of the eddies 
associated with the turbulent mixing of fuel and air. As stated previously, no data 
exist by which the significance of these non-uniformities and the role which turbulence 
plays on dissipating these non-uniformities has been generated for the case of autoig- 
nition of premixed systems. At best, theoretical models which treat the general class 
of turbulent chemical reactions are available; these models are reviewed in the follow- 
ing section of this report. 



7. Effect of Fuel Chemical Composition 

The chemical composition of a fuel can affect ignition delay time in at least 
two ways. First, the rate of fuel vaporization, and hence the fuel vapor concentra- 
tion, is determined by the volatility of the fuel; in the case of a multi-component 
fuel (e.g., a distillate fuel), the rate is determined by the relative proportions 
of high-volatility and low-volatility compounds in the blend. A procedure to cal- 
culate the vaporization rate of multi-component fuel blends is described by Anderson, 
et al. (1980). Second, the pre-ignition reaction system probably consists of a 
number of chemical reactions involving chemically simple molecules and radicals to- 
gether with one or more steps in which the large, complicated fuel molecules ther- 
mally decompose into these simpler structures. For example, Edelman, et al. (1972) 
proposed an ignition delay model based upon the assumption that the fuel decomposes 
into methane, ethane, and propane for which ignition delay reaction systems are 
known; this model was applied by Siminski and Wright (1972) to the analysis of the 
ignition delay data for a number of heavy hydrocarbon fuels. Hautman, ec al. (1981) 
present a model for the oxidation of paraffin-series hydrocarbons in which the hy- 
drocarbons are assumed to be converted to ethene which further reacts to form 
products of combustion. Empirical coefficients used in the model have been deter- 
mined using the results of shock tube measurements. A method of deducing the ig- 
nition delay period from the calculated behavior of the ether.e concentratione is 
proposed. Halstead, et al . (1975) have determined a rate limiting mechanism for 
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calculating the ignition delay time for a number of fuela used in reciprocating en- 
gines; they have determined the type of destruction reactions that limit the produc- 
tion of certain chain branching radicals. Affens and Carhart (1974) have attempted 
to relate the influence of chain branching radical distribution to fuel chemical 
characteristics. They have found that autoignition temperature increases with de- 
creasing number of carbon atoms in a hydrocarbon molecular chain and with increasing 
cyclicity or aromaticity; these findings indicate that ignition delay time increases 
with increasing complexity of the molecule's bond structure. 

It has also been determined by Ballal and Lefebvre (1979) that the minimum 
ignition energy for various liquid hydrocarbon fuels tends to increase with increas- 
ing molecular complexity and decreasing fuel volatility. It is possible to evaluate 
the relative importance of each effect from the data presented. 

Experimental studies have been carried out to determine the effect of fuel type 
and hence composition on ignition delay time. As indicated earlier, for those ex- 
periments conducted using liquid fuel sprays, it is not possible to separate fuel 
chemistry effects from physical effects (e.g., evaporation rate). Yashizawa, et al. 
(1978) conducted shock tube experiments using n-butane, n-hexane, and n-octane and 
concluded that, in the temperature range over which the experiments were conducted 
(900-1700K), there was no discernible difference in ignition delay time among these 
fuels. In comparing their results with other results, Yashizawa, et al. conclude 
that with the exception of the lighter hydrocarbons, methane and ethane, all satu- 
rated hydrocarbons exhibit nearly the same ignition delay values. Compared to the 
heavier saturated (paraffin-series) hydrocarbons, methane exhibits longer delay 
times while ethane exhibits shorter ignition delay times — times comparable to those 
observed for hydrogen and acetylene. These results suggest that differences in the 
ignition delay time behavior of liquid fuel sprays of the heavy hydrocarbon fuels 
will arise from the different physical characteristics alone; chemical effects can 
be ignored. 


ANALYTICAL METHODS FOR AUTOIGNITION MODELS 

On the basis of the experimental observations reviewed in the previous section, 
a model for the autoignition of fuel injected into a lean, premixing, prevaporizing 
fuel-air mixing passage must treat the occurrence of autoignition both in the vicin- 
ity of vaporizing droplets and in the fuel vapor-air mixture generated as the fuel 
vapor diffuses throughout the airstream. No analytic model was located during this 
literature survey in which this dual single-phase and two-phase autoignition problem 
has been modeled. However, a number of reports were obtained in which some impor- 
tant aspect of the overall problem has been treated. For convenience, the discussion 
of the results of the literature survey for analytical methods is divided into 
three parts: (1) a discussion of models in which no liquid phase is present, (2) 

a discussion of models applicable to the occurrence of autoignition in the droplet 
boundary layer, and (3) a brief discussion of turbulent reacting flow modeling. 
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Single Phase Models for Autoignition 

Two methods have been used to model the autoignition of fuel -air mixtures in 
which no droplets are present. In the first method, the reactants are assumed to be 
thoroughly mixed locally and a system of elementary reaction rates describing the 
pre-ignition process is assumed to apply. For example, it has been assumed that 
any large fuel molecule undergoes a one-step thermal decomposition process into 
various proportions of methane, ethane and propane [Edelman, et al. (1972)]; the 
calculation of ignition delay time is based upon a ten reaction rate system derived 
from autoignition data for these three simpler molecules. A similar model is des- 
cribed by Halstead, et al. (1975). However, in this case, the fuel molecule parti- 
cipates directly in the system of reactions by contributing hydrocarbon atoms and 
(upon partial oxidation of the fuel) hydroxyl radicals. Recently Hautman, et al. 
(1981) proposed a four-step kinetic mechanism for the oxidation of paraffin-series 
hydrocarbons in which the fuel first decomposes to ethene. The authors have shown 
that the onset of significant heat release corresponds to the time at which the 
ethene concentrations reaches a maximum; therefore a criterion serving to designate 
the ig.i ion delay period exits. All of these models have produced results in 
reasonable agreement with data for particular ranges of temperature and concentra- 
t ion . 


In the second approach to modeling autoignition in the absence of fuel droplets, 
the restriction that the fuel and air are locally premixed is removed [Dopazo and 
O'Brien (1974)]. Fuel and air are assumed to be contained in individual fluid 
eddies; intermixing of these eddies is controlled by the turbulence of the flow. 

The amount of contact time between the fuel and air is calculated from an assumed 
form of probability density function. The rate of reaction of fuel with air during 
the time intervals in which both reactants are mixed is estimated using a one— step, 
second order chemical reaction rate. Autoignition is indicated if the local tem- 
perature begins to increase rapidly, a phenomenon to which the authors refer as 
"thermal runaway". Thus, although the flow may appear to ’>e premixed on the basis 
of time-mean fuel and air concentrations, different reaction rates can be calculated 
depending upon the turbulence level and the assumed probability density function. 

Models for Autoignition in the Presence of Droplets 

Models for the autoignition of mixtures in which the fuel is injected initially 
as liquid droplets can be divided into two classes: (1) models that include the 

detailed structure of the mass and thermal boundary layers surrounding a droplet 
(referred to as in this report as boundary layer models), and (2) models that 
ignore these boundary layers except for the use of mass and heat transfer rates 
calculated a priori from knowledge of these structures. 

In the boundary layer models, the fuel vapor concentration and temperature pro- 
files in the boundary layer are calculated at each instant during the vaporization 
process. These profiles are used together with other criteria to determine the 
ignition delay time. In some cases, the pre-ignition period is divided into a number 
of subintervals with each subinterval characterized by a single rate controlling 
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process. For example, El Wakil and Abdou (1966) divided the pre-ignition period into 
physical delay and chemical delay subintervals. They defined the physical delay time 
as the time that elapses from the injection of the fuel droplet into the air stream 
to the development within the boundary layer of a fuel vapor concentration corres- 
ponding to the lean combustion limit. The chemical delay time is the time that elap- 
ses after the physical delay time until the occurrence of a rate of oxidation of fuel 
that is some fraction of the reaction rate at stoichiometric conditions. Calculation 
of the chemical delay time requires knowledge of a gas-phase global reaction rate for 
oxidation of the fuel. Satcunanathan (1971) presents a similar analysis which uses 
additional subintervals to define the ignition delay time. 

There are several boundary layer models in which the ignition delay time is not 
divided into subintervals. For example, Kadota , et al . (1976) first calculated the 
temperature and fuel vapor concentration profiles for the droplet and then estimated 
the ignition delay time at several positions within the boundary layers; the igni- 
tion delay time, x, is calculated as a function of the local temperature and fuel 
vapor concentration using correlations obtained from gas-phase experiments also 
described in their paper. The entire calculation is repeated for several times 
during the droplet vaporization process. Ignition is assumed to occur when 






o 

at some position within the boundary layers. This expression indicates that ignition 
has occurred when the residence time of the fuel-air mixture at a fixed position 
relative to the droplet center has exceeded the ignition delay time. With the addi- 
tional assumption that ignition is most likely to occur near stoichiometric fuel-air 
concentrations within the boundary layer, Law and Chung (1980) have developed a 
criterion for autoignition which specifies the relationship between a critical Dam- 
kohler number (rate of chemical reaction to rate of heat transfer) and parameters 
describing the droplet environment. 

In some boundary layer models, it is assumed that the pre-ignition reaction rates 
are controlled by the rate of production of an intermediate, but unknown, species 
X • for example, this spec es may be an important chain branching species, or it may 
be a species that destroys important chain branching species. Faeth and Olsen (1968) 
developed a model for autoignition in which the rate of production of the intermediate 
species is calculated using a single-step rate expression for the formation of Xi as 
a function of local fuel and air concentrations, temperature and pressure. The cal- 
culation is performed at various points within the boundary layer and several time 
intervals throughout the vaporization process. The rate expression is calibrated by 
adjusting the coefficients in the expression to give agreement with experimentally- 
derived correlations of ignition delay time data for the vapor of the same fuel in 
air. Ignition is assumed to occur when the normalized concentration of Xj (the con- 
centration of Xj divided by its value at the end of the pre-ignition period) exceeds 
unity. A similar approach has been used by Sangiovanni (1976) but the model described 
therein accounts for the diffusion of the intermediate species throughout the boundary 
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layer; under similar conditions, the effect of diffusion can be expected to increase 
the calculated value of ignition delay time relative to the value which is calculated 
using the model of Faeth and Olsen. It should be noted that since the models em- 
ployed by Faeth and Olsen and by Sangiovanni are calibrated using gas-phase ignition 
delay time correlations, they are essentially similar in approach to the method 
employed by Kadota, et al. (1976). 

There are autoignition models that do not consider the details of the fuel vapor 
concentration and thermal boundary layers surrounding the droplet. Instead, these 
models make certain simplifying assumptions and use expressions for the instantaneous 
droplet vaporization and heating rates; it should be noted that the derivation of 
these rates requires knowledge of the concentration and temperature profiles around 
the droplet. Rao and Lefebvre (1981) calculated the ignition delay time as the sum 
of a vaporization delay time and a chemical delay time. They defined the vaporiza- 
tion delay time as the time required to completely vaporize the fuel droplet; this 
time interval is estimated by the ratio of initial droplet mass to a constant value 
of calculated vaporization rate. Implicit in this model is the assumption that the 
pre-ignition reactions are unimportant until the fuel is completely vaporized. The 
rationale for the derivation of the expression for chemical delay time used in this 
model is ambiguous. Model constants are estimated from the available data. 

Turbulent Reacting Flow Modeling 

If the fuel vapor and air are mixed thoroughly on the molecular scale, then it 
should be possible to derive models for autoignition from an understanding of the 
important rate controlling processes and the available data; the assumption that 
reactions are thoroughly mixed is basic to the models presented by Edelman, et al. 
(1972) and by Halstead, et al, (1975) described earlier. If the flow is turbulent, 
then this assumption may not be justified. Each reactant may be contained in in- 
dividual eddies of fluid; the time-mean reactant concentration may indicate that the 
reactants are premixed, but the instantaneous concentrations may indicate that the 
reactants are unmixed for a significant portion of the time. If the latter situa- 
tion prevails, the calculated reaction rate will be less than that calculated for 
the premixed system. 

In modeling turbulent reacting flow, it is necessary to specify the dependent 
variables (velocities, temperature, pressure, species concentrations, etc.) for 
which fluctuations represent important departures from the mean values. For example, 
the instantaneous value of a dependent variable may be represented as u ® u + u'. 
Expressions for u together with similar expressions for the other dependent variables 
are then substituted into the equations of motion. After a number of algebraic 
manipulations that may include subtracting the corresponding equations for the mean 
values and neglecting certain higher order terms, one obtains a system of equations 
in terms of convective and diffusive terms for the fluctuating components together 
with products of two or more fluctuating components. These products are often 
termed correlations. Models for turbulent flow tend to be differentiated by (1) 
the number of fluctuating variables selected, (2) the order of the resulting equa- 
tions, and (3) the models used to calculate the correlations. A recent review of 
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turbulent flow models is presented by Libby end Williams (1981). Another review, * 
together with pepers on both enelyticel end experimental techniques in turbulent 
reacting flows is presented in e special issue of Combustion Science and Technology 
(1976). 

Most of the papers cited in modeling turbulent reacting flows le.g., Spalding 
(1971), Bilger (1976), Magnussen and Hjertager (1976), Lockwood (1977), Lockwood 
and Syed (1979) and Gosman and Ioannideo (1981)] are concerned with the rate of 
consumption of reactants. An exception is the model of Dopaso and O'Brien (1974) 
which is based upon solutions to the turbulent equations of motion using a statis~ 
tical approach for the modeling of the effects of the fluctuations in temperature 
of the reactants. During the pre— ignition period, negligible amounts of reactants 
are consumed, although the reactants that are consumed produce the intermediate 
products important to autoignition. It is difficult to extend these models for re- 
actant consumption to obtain an autoignition model because of the lack of data 
suitable for calibrating the model. For example, Magnussen and Hjertager propose 
that the rate of fuel consumption be calculated by an expression of the form: 



^ - AX f c/k (7) 

where e and k are determined from the solution of the turbulent equations of motion. \ 

The constant A varies from a value of 4 to 32 depending on the extent to which the 
reactants, fuel and air, are mixed. If it is assumed that the rate-controlling ’ \ 

intermediate species are produced at a similar rate (say) J 

* i 

(8) j 

j 

then data are required to determine B. Alternatively, a system of equations, analo- 
gous to Eq. (7), can be generated that are similar to the systems of equations des- 
cribing autoignition reaction rates in flows in which the effects of turbulence are 
neglected. The constants appearing in such a system would be adjusted to give J 

reasonable agreement with the available data. Unfortunately, no ignition delay 1 

time data which include the effects of turbulence have been published. I 



Apparently, attempts have been made to model the fuel oxidation process by 
starting with reaction rate expressions in the form (for example): 




into Eq. (9), turbulent forme of the reaction rate expreeaions are obtained. The 
fluctuating components of temperature and apecies concentration are obtained from 
the turbulent energy and species equations. Lockwood (1977) has stated that this 
approach is not useful because of the strong influence of temperature fluctuations 
due to the exponential term. This method can be useful in formulating a model for 
autoignition if (1) temperature fluctuations in the fuel preparation passage are 
negligible, (2) the appropriate pre-ignition reactions can be identified, and (3) 
data can be located relating ignition delay time to turbulence parameters. Since 
the system of rate equations must be solved at each point in the three-dimensional 
flow field, machine computation time requirements are likely to be so large as to 
make application of this approach impractical. 


CONCLUDING REMARKS 

The results of the literature survey indicate that the factors which are most 
significant in effecting the ignition delay time of fuel-air systems are the mixture 
pressure (reactant concentration), mixture temperature, and the presence of fuel 
droplets. Fuel composition and mixture ratio are of secondary impc tance. The 
role of turbulence has not been clearly defined. Analysis of the results of experi- 
ments in which liquid fuel is injected into hot flowing gas streams for the purpose 
of isolating the individual effects noted above is not possible because of the com- 
plexity of the flows and the resulting lack of knowledge concerning flow details 
(e.g., the initial spray droplet size distribution, the temperature profiles within 
the flow, etc.). The existing experimental results can best be utilized to support 
the development of the analytical model through comparisons of the analytical re- 
sults obtained when applying the model to the experimental conditions with the pub- 
lished experimental data. Obviously, a degree of engineering judgment will neces- 
sarily be employed in setting up the cases to be modeled because of lack of detail 
concerning the experiments; also judgments will be required to select which of the 
many empirical factors to be employed in the model should be changed in order to 
achieve agreement between experiment and theory. 


The results of this survey indicate that experimental data are lacking even for 
the most significant of factors controlling the process of autoignition of vaporizing 
fuel sprays. No definitive work on vaporization of fuel sprays has been performed. 
The kinetic mechanisms controlling autoignition of fuel at low (non-shock tube) 
temperatures is unknown. The role of the presence of droplet sprays (as opposed to 
single droplets or a series of individual droplets) is unknown. The analysis 
described in Section 8 reflects this stage of knowledge through the use of the 
most simple and * raight-forvard models of the individual processes as can 
reasonably be believed to account for observed experimental trends. 
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LIST OP SYMBOLS 

Conte *nts 

Droplet initial diameter 
Activation energy 
Reaction rate conatant 
Pressure exponent 
Pressure 

Universal gas constant 

Temperature 

Time 

Velocity, dependent variable 
Coneentrat ion 

Turbulent energy dissipation rate 
Turbulent kinetic energy 
Ignition delay time 

Fuel 

Intermediate species 

Referring to I, |, k chemical species 

Mean value 
Fluet unt Ion 
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FIG. A.2 


SUMMARY OF DATA CORRELATIONS SHOWING 
EFFECT OF TEMPERATURE ON IGNITION DELAY TIME 


F 

f 

t 

i 

t 

i 


l 

i 

\ 


a(t6.5) 


1000 


100 


10 


b(25 6) 


d(12 8) 



j(41 3) 


MAREK, ET AL (1977) 

SPADACCINI AND TeVELDE (1980) 
SPADACCINI (1977) 

SPADACCINI AND TeVELDE (1980) 
STRINGER. ET AL (1970) 

TABACK (1971) 

MESTRE AND DuCOURNEAU (1973) 
STRINGER. ET AL (1970) 

MULLINS (1953) 


FUEL 

JET-A 

NO 2 DIESEL 
JP-4 

NO 2 DIESEL 

AVTUR 

JP*4 

KEROSENE 

AVTAG 

KEROSENE 


J — i — I — i — I i I i I i l » I 


J L 


10 12 14 16 18 

1/T in , e ,x10 4 (K-1) 


20 


22 24 


L 


J L 




12 10 9 8 


Tjnlet x 10‘ 2 (K) 


• 1-11 - 1-3 


141 




